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ABSTRACT 


The  objective  of  this  thesis  v a s  to  examine  the  veather 
on  a  daily  basis,  with  a  viev  to  establishing  the  number  of 
vork  days  available  for  a  particular  farm  operation.  Although 
the  veather  may  influence  the  operation  of  farm  machinery  in 
several  "ays,  this  project  has  been  limited  to  consideration 
of  soil  moisture,  and  hence  soil  tractability  as  it  is 
influenced  by  the  veather. 

A  technique  developed  by  W.  Baier  and  G . ¥ .  Robert son 
vas  used  to  estimate  soil  moisture,  using  daily  minimum  and 
maximum  temperature,  and  precipitation,  obtained  on  IBM 
punched  cards  from  the  Meteorological  Branch,  Department  of 
Transport  in  Toronto.  In  addition  solar  radiation  at  the 
top  c f  the  atmosphere,  as  veil  as  monthly  averages  of  dev 
point,  v’ind  mileage  per  day,  hours  of  sunshine  per  day,  and 
d  a  y  1  e  rig  t  h  "ere  used. 

The  necessity  of  Including  the  above  mentioned  variables 
restricted  the  Alberta  stations  to  Medicine  Hat,  Lethbridge, 
Calgary,  Edmonton,  Vermilion,  Beaverlodge,  and  Fair view. 

At  some  stations  as  many  as  forty-five  years  vere  considered 
adding  up  to  a  total  of  102,000  station  days. 

In  order  to  facilitate  computer  handling  the  climatolog¬ 
ical  cards  vere  transferred  onto  magnetic  tape,  a  task  vhich 
involved  some  computer  Incompatibilities.  Once  on  tape  the 
data  vere  checked  for  mlspunched  data,  as  veil  as  for 
completeness  of  data,  and  a  complete  record  vas  produced. 
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Criteria  “based  on  traction  vers  developed  for  class¬ 
ifying  a  day  as  a  good  day  (vork  day),  or  a  bad  day  (non- 
vcrk  day).  Using  these  criteria,  and  the  amassed  veather 
data,  the  probability  of  a  bad  day  as  veil  as  the  prob¬ 
ability  of  a  good  day,  vere  determined  on  a  monthly  basis 
for  the  seven  Alberta  stations. 

It  vas  discovered  that  ’’persistence”  of  vork  days  or 
non-vork  days  does  occur.  This  means  that  the  probability 
of  a  vork  day  for  example  is  greater  if  the  previous  day 
vas  a  vork  day  than  if  the  previous  day  vas  a  non-vork  day. 
Hence  the  probability  of  a  day  being  a  vork  day  can  net  be 
established  objectively,  without  some  knowledge  of  conditions 
existing  on  the  previous  day. 

Although  the  probabilities  determined  in  this  project 
may,  or  may  not  be  entirely  reliable  for  all  combinations 
of  soil  and  topography,  they  do  help  to  quantify  variables 
that  are  of  great  significance  in  today’s  scientific  business 
of  farming.  Improvements  in  the  estimated  probabilities  for 
tillage  operations  as  veil  as  the  establishment  of  probability 


values  for  harvesting  operations  will  further  improve  the 
data  available  to  the  farmer  as  a  guide  to  correct  machinery 
sizing . 
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_ INTRODUCTION  AND  OBJECTIVES 

The  operation  of  farm  machinery  is  a  task  subject  to 
the  whims  of  nature,  and  optimum  machinery  sizing  requires 
some  consideration  of  weather  variability.  Before  attempting 
to  compensate  for  the  weather,  however,  first  consider  the 
manner  in  which  the  weather  may  halt  machinery  operations. 
Tillage  operations  are  In  general,  stopped  when  the  soil 
becomes  too  wet  to  provide  sufficient  traction  for  mobility, 
or  in  some  cases,  when  the  soil  becomes  vet  enough  to  stick 
to  the  implement.  Harvest  operations  may  be  stopped  in  an 
additional  manner,  in  that  the  soil  conditions  may  be 
satisfactory,  while  the  crop  may  be  too  vet  to  harvest. 

For  purpose s  of  this  thesis,  only  the  possibility  dealing 
with  traction  has  been  considered. 

The  objective  of  this  thesis  was  to  examine  the 
weather,  on  a  daily  basis,  with  a  view  to  establishing  the 
number  of  working  days  available  for  tillage  operations. 

The  general  problem  may  be  divided  into  several  sub- 
problems.  The  first  problem  was  that  of  establishing  the 
net  moisture  content  of  the  soil,  that  is  precipitation 
minus  evapo transpiration,  on  a  daily  basis.  Having 
established  the  moisture  content  of  the  soil,  it  was  then 
necessary  to  determine  the  relationship  between  soil 
moisture,  and  soil  tractability ,  in  order  to  classify  a 
particular  day  as  either  a  work  day,  or  a  non-work 


day.  When  this  relationship  was  satisfactorily  established, 
it  was  then  possible  to  determine  the  number  of  work  days 


. 
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available,  for  a  particular  tillage  operation.  The 
ultimate  goal  of  this  project  vas  to  establish  some  meaning 
ful  crobabilitv  values,  on  a  monthly  basis,  for  represent- 
ative  areas  of  the  province  of  Alberta. 


The  amount  of  available  literature,  dealing  vith  the 
general  problem  of  relating  cLirna ta logical  data  to  field 
tillability,  varies  depending  upon  the  particular  aspect  of 
the  problem  concerned .  Literature  dealing  vith  evaporation 
and  soil  moisture  content  is  relatively  abundant,  vhile 
literature  pertaining  to  the  relationship  between  soil 
moisture,  and  soil  trac tabili ty ,  or  implement  draft  require 
merits  is  relatively  scarce. 


2 . 1  Eva'po transpiration  and  Soil  Moisture 

Earlier  work  dealing  with  evaporation  was  concerned 

15 

vith  rater  loss  from  open  water  surfaces.  Chew  reported 

chat, ’’Under  given  conditions,  evaporation,  is  proportional 

to  the  deficit  in  vapor  pressure,  which  is  the  difference 

between  the  pressure  of  saturated  vapor  at  the  tempera ture 

of  the  water  and  the  aqueous  vapor  pressure  of  the  air." 

Whether  over  open  water,  or  over  soil,  evaporation  is 

related  to  the  energy  balance  at  the  earth’s  surface. 

19 

Halstead  and  Covey  note  that  this  balance  depends  upon 
four  factors:  the  radiation  received,  the  heat  flux  to  the 

soil  by  conduction,  the  heat  flux  to  the  air  by  convection, 
and  the  heat  of  vaporization  of  water.  The  apportioning 
of  net  radiation  into  the  remaining  components  depends  upon 
a)  the  temperature  distribution  in  the  soil  and  in 
the  air  and  the  water  vapor  distribution  in  the 


air . 
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and  on  the  rate  of  movement  and  mixing  cf  air  in 
and  over  the  plant  cover. 

Energy  relat ionships ,  and  soil  moisture  retention, 
involve  the  concept  of  capillary  potential,  which  is  the 
amount  of  work  required  to  pull  a  unit  mass  of  water  from 


a  unit  mass  of  soil.  This  idea  is  explained  with  the  use 

15,  6  ‘  55 

cf  calculus  by  Chow  ,  Baver  ,  and  Richards  .  Richards 

56 

et  al.  simplified  this  theory,  and  expressed  soil 

moisture  content  as  a  function  of  time  after  irrigation. 

Because  of  the  complexity  of  the  energy  balance  approach, 

simpler  means  of  estimating  evaooration  were  sought. 

42 

Thornthwaite  (1948)  pioneered  this  line  of  thought  by 
estimating  evaporation  from  water  surfaces,  using  the  mean 
monthly  temperature,  and  a  heat  index  which  varied  from 


31 

station  to  station.  Penman  (1948)  developed  a  more 
sophisticated  means  of  estimating  evaporation  using  wind 
velocity,  vapor  pressure,  solar  energy,  and  temperature. 
This  method  it  more  accurate  than  Thornthvai te 1 s . 

The  closeness  of  the  relationship,  between  evapo tran¬ 
spiration  and  solar  radiation,  and  the  necessity  of 

including  radiation  in  estimates  of  evanoration  was  dis- 

44  30 

cussed  by  Pruit  and  Van  Vi jk  .  Pel ton  reported  that 

between  70  and  80$  of  the  evaporation  from  BellanJ.  plates, 


Class  A  veather  pans,  and  the  experimental  farm's  buried 
tank  at  Swift  Current,  could  be  explained  by  air  temp¬ 
erature,  radiation,  wind,  and  vapor  oressure.  At  Summer land , 

45 

B.C.,  Wilcox  found  a  high  correlation  between  evapc tran¬ 
spiration  from:  alfalfa,  ladino  clover,  and  brome  grass 

and  solar  radiation,  temperature,  wind,  and  humidity.  In 

20 


Southern  Alberta,  Hobbs  and  Krogman  correlated  evap- 
c transpiration  estimates  by  Thcrnthvsite  and  Penman  with 

evapo transpiration  from  second  year  alfalfa.  Hobbs  and 

21 

Krogman  also  related  evapo transpiration  from  alfalfa  to 

meteorological  variables.  Thus  the  soundness  of  using 

clima ta logical  factors  to  estimate  evaporation,  and  potential 

evapo transoiration,  has  been  well  established. 

25 

Lignon  et  al .  used  Thornthvaite * s  method  of  estimating 

potential  evapo transpiration,  which  they  combined  with 

precipitation  to  determine  soil  moisture.  To  improve  their 

estimates  they  added  corrections  for  soil  dryness  and  for 

the  Influence  of  the  veather  on  evapo transpiration.  A 

52 

similar  method  was  used,  by  Pierce  ,  -ho  used  mean  monthly 
temperature,  with  corrections  for  soil  dryness,  crop  stage, 
rainy  days,  as  rell  as  day  length. 


5 

Baier  and  Robertson  used  maximum  temperature,  temp¬ 
erature  range,  solar  radiation  at  the  top  of  the  atmosphere, 
hours  of  sunshine,  day  length,  wind  mileage,  and  vapor 
pressure  to  estimate  the  daily  latent  evaporation  from 


Bellani  plates.  They  obtained  a'  correlation  coefficient 
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a  no 
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22  24 

of  0.84.  Ho  lines  and  Robertson  ,  Krogman  and  Hobbs, 

Baier  converted  latent  evaporation  to  potential  evap- 
o transpiration  by  multiplying  latent  evaporation  by  0.00^4. 
Using  potential  evapo transpiration  .  and  precipitation, . vith 

corrections  for  crop  stage  and  soil  dryness,  Baier  and 

4 

Robertson  .estimated  soil  moisture  under  vheat  for  various 

6 

soil  depths.  Baier  and  Robertson  found  a  high  correlation 
betveen  estimated  soil  moisture  and  crop  yields  for 


several  points  throughout  Canada. 

The  idea  of  estimating  soil  moisture  using  a  regression 
of  clima talogical  variables  is  seen  to  be  veil  supported. 
Because  Baier  and  Robertson's  technique  appeared  to  be  the 
most  refined  it  was  used  here. 

2 .2  Soil  Moisture  and  Sol 1  Tra c tability 

The  specific  problem  of  relating  soil  tractability  to 

soil  moisture,  vith  special  regard  to  farm  machinery,  does 

not  appear  to  have  been  examined.  The  problem  of  tractability, 

and  soil  bearing  capacity  under  veak  soils  has  been  examined 

.in  connection  vith  the  design  of  earth  crossing  vehicles 

9,  10  4p 


by  Bekker  and  by  the  U.S.  Army  .  The  general  problem 

of  soil  strength  and  soil  tractability  has  been  examined 

9,  10  23  46 

by  several  vriters  (Bekker  ,  Janosi  ,  Wills  ,  Barger 

7  34 


et  al .  ,  Reece  ).  These  vriters  expressed 
resiscance,  and.  soil  thrust  in  terms  of  the 
variables:  cohesion  of  soil  (c),  internal 


s i nk age,  rolling 
fol loving  soil 
frictional  angle 
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of  noil  W-) ,  cohesion  modulus  of  soil  deformation  (Kc)  , 
frictional  modulus  of  soil  deformation  (Xf/),  and  the  coeffic¬ 
ient  of  slnkage  (n),  These  soil  variables,  must  be  determined 
experimentally,  for  soils  of  varying  texture,  and  moisture 
contents  before  the  relationships  described  by  these 
writers  can  be  used  in  relating  soil  moisture  to  soil 

tree  Lability , 

'* 

The  shear  strength  of  a  soil  is  given  by  Coulomb’s 
equation,  S  -  c  p  tan  where  p  is  the  confining  pressure 

and  c  and.  are  as  defined  above.  The  angle  of  internal 

17 

friction  in  top  soils  vas  found  by  Founts ine  and  Brown 

to  be  quite  high  under  small  normal  loads  (o-lO  psi). 

16 

Evans  .  indicates  that  the  shear  strength  of  dry  top  soils 

depends  upon  their  compaction.  However,  if  the  moisture 

content  is  well  above  the  plastic  limit,  compaction  reduces 

the  voids  causing  the  soils  to  become  saturated.  In  this 

case  compaction  does  not  increase  soil  strength.  He 

found  that  agricultural  loams,  with  a  clay  content  above 

20%,  and  a  moisture,  content  in  the  plastic  range,  appear 

to  have  no  angle  of  shear  resistance.  Sandy  soils  however, 

have  a  high  angle  of  shear  resistance,  and  no  cohesive 

37 

strength.  Skempton  and  Bishop  indicate  that  both  cohesion 
and  interna l  friction  depend  upon  the  void  ratio  of  the 
soil  at  the  instant  of  failure.  Hence,  the  determina t ion 
of  soil  trac tab ility  is  seen  to  be  quite  complicated. 


-8 


29 

Nichols  experimented  vith  four  clays,  a  clay  loam, 
and  a  silt  loan,  and  established  a  relationship  be  m  ean 
clay  content,  liquid  limit,  moisture  content,  confining 
or e s s ure ,  a  nd  sh a a r  s tr e ng th  of  pi a s tic  s o i 1 s  .  S in c e  the 
values  required  using  this  technique  -ere  available  for  a 
number  of  soils  in  the  province,  this  "as  the  technique 
used  t  o  esc i ma t e  s oil  shear  s tr e ng  th  i n  th i s  s  tud y . 

In  addition  to  kncving  the  shear  strength  of  the  soil, 
an  estimate  of  the  shear  strength  necessary  to  provide 
traction  vas  also  required,  -  This  is  determined  by  the 
draft  required  to  pull  a  tillage  implement,  the  rolling 
resistance  of  the  tractor,  and  the  contact  area  of  the 
tires  vith  the  soil,  he  literature  mas  found  relating  any 
of  these  factors  to  soil  moisture  content.  Hovever,  rolling 


resistance  has.  been  related  to  outside  tire  diameter 


7,  ^ 


Draft  requirements  for  various  Implements  under  general 

r7 


conditions  are  given  by 


.r.v,g 


r  e  t  ai .  a  nd  i n  A . M . A .  t  e  s t 


A 

reports  . 


While  information  dealing  vith  specific  topics  was  In 
some  instances  not  precise,  some  information  at  least  vis 
found  dealing  vith  the  topic  of  soil  tracts bility  in  general 


* 


SOIL  MOISTURE 


2_. 

Net  soil  moisture  change  depends  upon  two  quanties: 
precipitation  and  evapo transpira tion.  This  moisture  change 
was  estimated  using  a  technioue  developed  by  Baler  and 

"4  - 

Robertson  .  A  resume  of  their  technique  supplemented  by 
private  communications  as  applied  to  Alberta  conditions 
is  in  order. 


3 . 1  The  Versatile  Soil  Moisture  Budget 

In  the  Versatile  Soil  Moisture  Budget,  the  soil  was 
divided  into  six  soil  moisture  zones,  holding  0.20  in., 

0.30  in.,  0.50  in.,  1.0  in.,  1.0  in.,  and  1,0  in.  of  water, 
respectively.  These  zones  were  based  on  a  fixed  depth  of 
water  held  in  the  soil,  rather  than  on  a  fixed  depth  of 


soil,  in  an  attempt  to  make  the  budget  independent  of  soil 
type.  Hence,  the  thickness  of  the  soil  layer  concerned, 
depended  upon  its  moisture  holding  capacity. 

The  addition  of  moisture  to  the  soil  was  for  the  most 
part  equal  to  the  precipitation.  An  exception  to  this  genera 
case  was  a  twenty-four  hour  rainfall  in  excess  of  one  inch. 

In  this  case  it  was  assumed,  that  some  run  off  occurred. 

If  run  off  did  occur,  it  was  equal  to  precipitation  minus 
infiltration  where  infiltration  (l)  is. 


1 


I  ^  0.9177  +  1.811  loge  (rainfall)  T  0.0097  loge 
S  ‘j  (i-lji 

(rainfall)  s.j  .  100,  where  S.j  is  the  capacity 

of  zone  1  and  S'j(i-l)  is  the  available  moisture 
at  the  end  of  the  day  preceding  the  rainfall. 


■ 
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On  rainy  days,  -evaporation  vas.  assumed  to  have  occurred 
in  the  morning,  while  precipitation  was  assumed  to  have 
occurred  in  the  afternoon.  If  the  rainfall  was  associated 
with  a  low  pressure  centre,  or  due  to  orographic  lifting 
as  a  result  of  a  high  pressure  centre,  this  assumption  may 


he  invalid,  hut  as  is  the  usual  case  in  Alberta,  if  the 
rainfall  was  of  a  showery  nature  the  assumption  is  reasonable 

In  the  estimate  of  evapo transpiration,  allowance  was 
made  for  stage  of  crop  development .  This  was  accomplished 
by  the  use  of  the  coefficients  given  in  Table  1. 

The  coefficients  for  the  period  from  planting  to  the 
hard  dough  stage  were  the  best  available  for  estimating 
soil  moisture  under  wheat.  Because  of  a  desire  to  extend 
the  budget  beyond  this  range,  the  coefficients  for  the 
period  from  April  1  to  planting,  and  the  period  from  hard 


dough  to  October  31*  were  estimated  by  the  author.  The 
coefficients  for  the  period  from  April  1  to  planting  were 


assumed  to  be  the  same  as  from  planting  to  emergence, 
whereas  the  coefficients  from,  hard  dough,  to  October  31,  were 
assumed  to  be  between  those  for  the  period  planting  to  emergence, 
and  the  period,  soft  dough  stage  to  hard  dough  stage.  The 
latter  coefficients  can  no  doubt  be  improved.  The  develop¬ 
ment  of  similar  coefficients,  describing  the  uptake  of  water 
by  other  crops,  as  veil  as  the  loss  of  water  by  a  summer - 
fallow  field,  would  also  improve  the  usefulness  of  the 
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technicue,  but  for  present  purposes,  only  moisture  conditions 
under  wheat  could  be  considered. 


Tabl 
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ef 
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on  the  relationship  between  potential  evapo transpiration, 

and  actual  evapo transpira tion.  The  nature  of  this  effect 

may  vary  under  different  environmental  conditions.  Figure 

1  indicates  various  AE/PE  relationships  as  proposed  bv  Baier 

4 

and  Robertson  .  One  or  more  of  the  proposed  AS/PE  relationships 


could  he?. ve  been  incorporated  into  the  program  for  estimating 
soil  moisture  by  including  the  required  coefficients  in  the 
form  of  ten  by  ten  matrices.  Tables  2  and  3  include  the 


Coefficients  for  planting  to  hard  dough  were  obtained  in 
personal  communication  vith  G.  W. 


Robertson . 
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coefficients  describing  curves 
of  relationship  applicable  to  v 
at  Ottawa.  The  use  of  curve  D 


D  and  E.  Curve  E  shows  the  type 
heat  grown  .In  Mathilda  loam 


after  harvest  at  Otta 


M 


cuced 


the  discrepancies  between  observed  and  estimated  soil  moisture 
for  the  0.20  inch  zone,  having  no  basis  for  doing  otherwise, 
curve  S  ’-'as  usee?  before  harvest,  while  curve  D  was  used 
after  harvest.  Perhaps  different  soils,  in  different  areas 
of  the  province  require  somewhat  different  tables  of 
coefficients,  but  because  there  is  no  experimental  work 
for  the  province  the  curves  used  by  Baler  and  Robertson 
were  used  here . 


- 
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Table  2, *  Ad ,j us t me n  t  Factors ( Z )  for  Soil  Dryness  Curve  D 
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In  the  Versatile  Soil  Moisture  Budget,  water  was  with¬ 
drawn  simultaneously  from  different  depths  of  the  soil 
profile  in  the  root  zone  in  relation  to  the  potential 

*  Coefficients  vers  obtained  in.  personal  communication 
v i th  G .  W .  Hob er t  so  n 


• 

• 
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cvaootransolra tion  rate  and  the  available  soil  moisture  in 

<L 

each  zone.  The  actual  daily  evapo'transpira  t  Lon  vs  s  estimated 
from  the  following  equation: 


a  .a  l  - 


where 


l 

.1  =  1 


v-  -  Q| 


k.) 


s" 


.1  ( i  ~  1)  .  Zj.PEi*9-v(PEi  r  FI?) 
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summation  carried  p.ut  from  zone,  j  »  1  to.  zone  j  »  n 
~  coefficient  for  soil  and  olant  characteristics  in 


zone  j 


q  i  if 


j(i  -  1) 


available  soil  moisture  in  the  jth  zone  at 
the  end  of  day  i  -  1 


/ 1  * 
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PS  rates  on  AS/pE  ratio 
F?  =  average  monthly  PF 


Potential  evaoo transpiration  has  been  estimated  by 
3,  22,  24 

several  corkers  as  being  equal  to  the  latent 

evaporation  (Li )  times  0.0034.  Latent  evaporation  is 
either  measured  with  a  Bellani  plate  atmometer,  or  it  may 
be  estimated  from  weather  parameters,  as  was  the  case  here. 
The  technique  used  in  estimating  latent  evaporation  is 
discussed  in  section  3.2 
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The  term  e~v  (i>rV^  “  accounts  for  varying  atmospheric 

demand  rates  (pSi)  on  the  relative  evapo transpiration,  as 

compared  to  the  long  term  average  PE  for  the  month  (PE),  under 

different  soil  moisture  stresses  (v ) .  Soil  moisture  stress 

is  dependent  upon  the  available  soil  moisture  in  a  particular 

zone,  and  is  expressed  as  a  percent  of  the  available  water 

holding  capacity  of  that  zone.  A  regression  equation  (r  = 

0.866)  used  for  estimating  w,  has  been  developed  by  Baier 

4 

a nd  Rober t so n  ,  as  follows: 

v  =  7.81  -  0.11  .  SM(i  -  1).  100. 

Sj 

The  soil  Moisture  Budget  required  some  initial  conditions 
and  controls,  whether  they  be  observed  or  estimated.  In 
this  case  they  had  to  be  estimated. 

Monthly  average  potential  evapotranspiration  rates  were 
required  for  the  months  of  concern  for  each  station  con¬ 
sidered.  These  rates  were  determined  by  the  technique 
discussed  in  section  3.8  and  are  found  in  Table  4. 

An  estimate  of  soil  moisture  at  the  beginning  of  each 
season  was  also  required.  The  assumption  made' was  that 
the  moisture  content  of  each  zone  was  at  one  half  capacity 
at  the  beginning  of  the  first  season  for  each  station.  The 
moisture  content  at  the  beginning  of  each  successive  season 
was  estimated  by  adding  a  fraction  of  the  winter  precipitation 
to  the  fall  moisture  content.  The  fraction  of  winter 

precipitation  that  is  absorbed  by  the  soil  varies  with 

4l 

latitude.  Staple  and  Lehane  of  Swift  Current,  found  that 
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between  November  and  April,  between  1940  -  1950,  330  of  the 

14 

precipitation  was  conserved.  Dr.  Carder  at  Beaverlodge 

found  that  from  harvest  to  freeze  up  130  of  the  precipitation 

was  conserved  whereas  from  freeze  up  till  spring  break  up 

50  was  conserved.  For  purposes  of  this  study,  it  was  assumed 

that  for  stations  (the  choice  of  stations  is  discussed  in 

section  3.5)  Medicine  Hat,  Lethbridge,  and  Calgary,  3 5$  of  the 

winter  precipitation  was  conserved,  at  Edmonton  and  Vermilion 

25^  vas  conserved,  and  at  Beaverlodge  and  Fairviev  100 

was • conserved .  The  amount  of  precipitation  conserved  in  the 

southern  part  of  the  province  may  be  far  less  than  at  Swift 

Current  because  of  the  influence  of  Chinooks,  but  no 

information  was  found  for  Southern  Alberta  stations. 

Crop  stage  dates  were  also  required.  Upon  examination 

o 

of  stage  dates  (Baier  and  Robertson  )  for  Beaverlodge, 


La combe,  and  Swift  Current,  there  appeared  to  be  little 
difference  in  stage  dates  from  one  station  to  another,  until 
the  soft  dough  stage,  after  which  a  difference  of  two  weeks 
was.  found.  As  a  matter  of  convenience  stage  dates  were 
assumed  to  be  constant  throughout  the  province  and  they  were: 
emergence  May  20,  jointing  June  17,  heading  July  14, 
flowering  July  19,  milk  August  2,  soft  dough  August  10, 


and  hard  dough  September  7. 

The  usefulness  of  the  program  would  be  increased 
markedly  if  provision  were  made  for  variation  in  planting 


* 
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dates,  and  consequent  later  stage  dates.  The  stage  dates 
following  planting  are  dependent  upon  a  number  of  factors. 
Including  precipitation  and  hours  of  sunshine.  Because 
these  dates  are  not  easily  determined  average  dates  had  to 
be  assumed  for  these  stages.  If  allowance  could  be  made  for 
these  variations  the  program  would  obviously  be  of  much 
gr e a  t e r  bene fit. 


2^.2 _ Estimation  of  Latent  Evaporation 

The  Soil  Moisture  Budget  required  an  estimate  of  latent 

evaporation,  which  vas  determined  separately  by  Baier  and 
5 

Robertson  .  Using  three  parameters,  minimum  and  maximum 
temperature,  and  solar  radiation  at  the  top  of  the  atmosphere, 
they  found  a  correlation  coefficient  of  0.58.  By  including 
additional  factors  of  sunshine  hours,  day  length,  dew 
point,  and  wind  mileage,  the  correlation  coefficient  vas 
raised  to  0.84.  In  other  words,  JQ$>  of  the  variation  was 
explained  by  variation  in  the  above  factors. 

Either  daily  observations  or  long  term  monthly  averages 
of  wind  mileage,  hours  of  sunshine,  dev  point,  and  day 
length  could  have  been  used.  Although  not  as  accurate  as 
daily  observations,  monthly  averages  to  help  make 
allowance  for  station  weather  differences.  Latent  evaporation 
vas  estimated  from  the  following  equation, 

EL  -  -53.39  +  0.337MAX+  0.531  (MAX  -  MIN)  + 

0.0512QS  +  0.0977  WIND  +  1.77.  (ev  -  es) 


0.0107QO 


• 
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vhere 

El  -  latent  evaporation 

MAX  -  maximum  Fahrenheit  temperature 

MIN  -  minimum  Fahrenheit  temperature 
Qo  =  solar  radiation  at  the  top  of  the  atmosphere 
Qs  -  solar  radiation  received  on  a  horizontal  surface 
WIND  =  total  daily  wind  mileage 

ey  -  saturation  vapor  pressure  at  mean  air  temperature 
e3  -  saturation  vapor  pressure  at  mean  dev  point 

The  term  EL  indicates  evaporation  from  a  Be-llani  plate 
atmometer.  It  is  measured  in  cc.  per  day  and  in  this  case 
was  an  estimated  value. 

Minimum  and  maximum  temperature  were  both  recorded  in  a 
Stevenson1 s  screen  four  feet  above  ground  level.  These 
values  were  obtained  from  Type  4  Climatalogical  cards,  which 
were  obtained  from  the  Meteorological  Branch,  Department  of 
Transport,  in  Toronto. 

The  entity  Qo  is  the  solar  radiation  at  the  top  of  the 

atmosphere.  Values  of  Qo  were  obtained  on  IBM  cards,  one 

card  per  day  for  each  degree  of  latitutde,  from  the  Plant 

Research  Institute,  Ottawa.  They  are  also  available  in 

38 

Smithsonian  Tables 

The  term  Qs  is  the  solar  radiation  received  at  the 
earth.  It  was  derived  from  a  separate  regression  equation, 
with  a  coefficient  of  determine t ion  of  0,85.  The  equation 
also  developed  at  Ottawa  is  as  follows: 
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Qs  =  Qo( 0.251  +  0.6l6(n/N)) 

where  n  is  the  daily  hours  of  bright  sunshine  and  N  is  the 
number  of  hours  between  sunrise  and  sunset.  The  value  n 
was  taken  as  a  monthly 'average  (the  reason  being  explained 
in  section  3*3)  and  a  summary  of  n,  for  each  month  at  each 
station,  is  found  in  Table  5.  In  this  case  a  monthly  average 
was  used.  The  monthly  averages  were  obtained  from  Smithsonian 
tables  and  are  summarized  in  Table  6. 

Wind  mileage  was  also  a  monthly  average.  The  anemometer 
at  Medicine  Hat  was  located  57  ft.  above  ground,  at  Lethbridge 
85  ft.  (changed  to  10  m.  November,  1965),  at  Calgary  57 
ft.,  at  Edmonton  67  ft.,  at  Vermilion  38  ft.,  at  Beaverlodge 
35  It.,  and  at  Fairview  38  ft,  above  ground.  Because  of  the 
difference  in  anemometer  heights,  a  wind  reduction  formula 


vas  us 

ed 

to 

bring  wind  mi 

leage 

formu'l 

a  .1 

log  h£ 

u6 

= 

ux  . 

log  hx 

where 

Ux 

S 

wind 

at  height 

hx 

% 

wind 

at  6  feet 

-o* 


This  formula  was  suggested  by  Baier  and  Robertson  (private 
communication) .  Summaries  of  the  wind  mileage  for  the 
different  stations  used  are  given  in  Table  7, 

Saturation  vapor  pressures  were  obtained  from  tables 
The  dev  point  vas  again  a  monthly  average.  The  dev  point 
for  each  station  is  given  in  Table  8, 


38 
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3.3  Programing 


'The  determina tion  of  latent  evaporation,  and  sub¬ 
sequently  soil  moisture  is  a  lengthy  process,  however 
computations  are  greatly  facilitated  by  the  use  of  a  digital 
computer.  At  the  time  this  project  was  started  computer 
programs  written  in  Fortran  II  were  obtained  from  the 
Plant  Research  Station  in  Ottawa  'go  compute  both  latent 


evaporation  and  soil  moisture.  Programs  are  now  available 
in  Fortran  IV  from  Ottawa. 

The  program  for  estimating  daily  latent  evaporation 
would  accept  either  daily  observations  or  monthly  averages 


of  wind  mileage,  dev  point,  hours  of  sunshine,  and  day 
length.  In  addition,  daily  observations  of  temperature  and 
precipitation  were  required  as  input.  The  output  from  this 
program  either  in  the  form  of  punched  cards  or  on  magnetic 
tape  contained  the  date,  estimated  latent  evaporation,  and 
precipitation.  This  information  then  provided  the  input 
for  the  soil  moisture  program.  Both  of  these  programs  were 
much  more  versatile  than  vas  required,  and  since  both  had 
to  be  rewritten  in  Fortran  IV,  the  two  programs  were  greatly 
simplified  and  combined  into  one  program. 

Since  monthly  averages  were  used  for  all  the  weather 
parameters  except  temperature  and  precipitation,  the  combined 
program  has  provisions  for  accepting  daily  observations  of 
temperature  and  precipitation  only.  Hence  the  versatility 
of  the  programs  vas  reduced,  but  soil  moisture  vas 


calcula ted 


- 


. 
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f  rom 


stop 


the 


of 


basic  weather  parameters  without 
obtaining  evaporation  output  and  f 


the  in  t errne d  i a  t e 
eeding  it  into  a 


s  o i L  moist ur e  p  r o gr a  m . 

The  soil  moisture  content  T-~as  used  to  determine 
whether  tillage  operations  were  possible  or  not.  A  good  day 


was  represented  by  an  X  and  a  bad  day  was  represented  by 
a  dot  (.),  with  one  line  of  X’s  and  dots  representing 
a  time  interval  of  one  month.  Hence,  the  combined  program 
accepts  as  input  Type  4  Climatalogical  cards,  and  summarize 
the  output  in  the  form  of  one  line  per  month,  indicating 
field  workability.  If  more  versatility  in  the  soil  moistur 
program  is  required  the  revised  programs  from  Ottawa  must 
be  referred  to. 


3 . 4  _ Selecting  and  Obtaining  Data 

The  techniques  developed ‘by  Baier  and  Robertson  required 
daily  observations  of  minimum  and  maximum  temperature,  and 
precipitation.  Also  required  were:  observations  of  wind 
mileage,  hours  of  sunshine,  and  dew  point,  either  daily, 
or  as  monthly  averages.  Hence,  the  problem  determining  the 
data  required  was  solved. 

Upon  contacting  the  Meteorological  Branch,  Department 
cf  Transport,  Toronto  it  was  learned  that  the  required 
information  was  available  on  IBM  punched  cards.  However, 
the  desired  daily  data  were  not  all  readily  available  on 
single  cards.  Instead,  minimum  and  maximum  temperature,  and 
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precipitation,  vere  available  on  one  card ;  vind  mileage  vas 
available  on  a  second  card;  hours  of  sunshine  vas  available 
on  a  third  card;  and  de*  point  vas  available  on  yet  another 
card.  Wind  velocity  and  dev  point  data  ver  punched  on  an 
hourly  basis,  and  while  cards  could  have  been  obtained  .for 
selected  hours,  the  cost  vould  have  been  as  great  or  greater 
than  if  the  cards  for  all  tventy-four  hours  vere  taken. 
Irrespective  of  cost,  if  fifty  cards  per  day  (the  number 
required  to  obtain  all  the  required  information)  had  been 
obtained  for  a  number  of  years  at  a  number  of  stations, 
the  volume  of  cards  required  would  have  been  astronomical. 
Another  costly  alternative  vould  have  been  to  have  obtained 
a  composite  card  containing  all  the  desired  information. 

.Jot  only  vould  this  have  been  costly,  but  a  ’waiting  period 
of  .from  five  to  six  months  vould  have  been  required.  It 
vould  also  have  been  possible  to  have  obtained  the  data  on 
tape . 

Partly  because  of  cost  considerations,  and  partly 
because  the  moisture  budgeting  technique  of  Baler  and 
Hob  .rtson  vould  accept  monthly  averages,  only  the  cards 
co nt a  ini ng  ;n i n i mum  a nd  rna x  1  m urn  temperature,  and  p r  ' c  L p i t a  tion 
(Type  4  Climatological  Cards)  ’.  ere  obtained.  These  cards 
were  available  at  a  cosc  of  six  dollars  p  r  thousand  cards. 
The  remaining  data  ’were  obtained  from  monthly  summaries. 

The  monthly  averages  of  vind  mileage,  hours  of  sun¬ 
shine,  dev  point,  and  potential  evapo transpiration  summarised 


~x. 

>  — 


in  Tables  4  - 


U 


ere  based  on  a  ten  year  period  1957  - 


i960,  for  the  months  April  to  October.  For  the  remaining 
months,  the  values  for  1966  vers  used  simply  to  provide  a 
value  in  order  to  make  the  program  vork,  since  evaporation 
vas  of  little  interest  during  these  months. 


5 .5  The  Selection  of  Stations 

The  choice  of  stations  from  vhich  to  obtain  veather 
data  vas  based  primarily  upon  the  number  of  veather  par¬ 
ameters  recorded  at  individual  stations.  Since  observations, 
or  monthly  averages  of,  minimum  and  maximum  temperature, 
precipitation,  v ind  mileage,  hours  of  sunshine,  and  dev 
point  vere  required  for  estimations  of  evaporation,  and 
subsequently  soil  moisture  content,  the  stations  had  to  be 
chosen  from  the  stations  that  recorded  all  these  parameters. 
Medicine  Hat,  Lethbridge,  Calgary,  Edmonton,  and  Fairviev 
vere  the  only  stations  in  .the  province  that  recorded  all  the 


required  parameters  for  a  sufficiently  long  period  of  time. 
Beaverlodge  recorded  all  the  parameters  excepting  dev 
point.  Since  Grande  Prairie  recorded  dev  point  and  vas 
located  only  thirty  miles  from  Beaverlodge,  it  as  felt 
that  the  dev  point  at  Grande  Prairie  vould  be  2  sufficiently 
close  estimate  of  dev  point  at  Beaverlodge.  Vermilion 
did  not  record  hours  of  sunshine,  but  Ranfurly  did,  and 
again  it  vas  felt  that  conditions  at  Vermilion  vould  be 
sufficiently  close  to  conditions  at  Ranfurly,  to  permit  the 


Table  4.  Monthly  Average  potential  S  va  po  tr  a  asp  Ira  t  ic-  n 
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Table  7  Monthly  Average  Dally  Wind  Mileage  (at  six  feet 
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use  cf  sunshine  at  Ranfurly.  With  these  assumptions  in 
mint?,  composite  stations  vere  formed  at  these  c*o  stations. 

The  fact  that,  the  number  of  stations  available  vas 
limited  made  it  impossible  to  select  stations  at  random. 

In  fact  every  suitable  station  vas  chosen.  If  the  province 
were  divided  up  into  its  component  farming  areas,  samples 
representing  both  farming  areas  and  geographic  areas  of 
the  province  could,  be  obtained  by  choosing  a  station  from 
each  area . 

Figure  2  can  be  used  to  point  out  the  various  farming 
areas  of  the  province.  The  map  -as  prepared  by  V.  Anderson  . 

The  drier  ranching  area  of  the  province,  designated 
as  region  A,  has  very  little  crop  production.  This  area 
of  the  province  vas  justifiably  excluded  from  a  study 
involving  the  influence  of  veather  on  machinery  operation, 
and  has  been  emitted. 

Region  B  which  includes  most  of  the  irrigation  area 
of  Alberta,  plays  a  large  economic  role  in  the  farming 
segment  of  the  economy.  Although  irrigation  farming  has 
not  been  considered,  dry  land  farming  in  that  area  of  the 
province  vas  represented  by  Lethbridge. 

A  large  grain  producing  area  of  the  province,  designated 
as  region  C,  “as  represented  on  opposite  corners  by  Calgary, 
and  by  Medicine  Hat.  Except  along  the  western  edge,  yields 
are  not  high,  but  large  acreages  compensate  for  lov  yields, 
hence  it  was  fitting  that  this  area  be  represented. 
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Figure  2.  Map  of  Alberta  ShoT.  ing  the  Division  of  Alberta 
into  the  Predominant  Farming  Patterns  and  The 
Location  of  the  Weather  Stations. 
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Region  D  is 


the  foothills  ranching  area 


of  the  province. 


Being  unsuited  to  crop-  production,  its  exclusion  from  the 
study  was  justified. 

Edmonton  and  Vermilion  each  represent  a  mixed  farming 


area  of  the  province  labelled  region  E.  Acreage  per  farm 
is  less  than  in  region  C,  but  capital  investment  per  farm 
is  comparable.  Because  a  large  percentage  of  the  land  in 
this  area  is  cultivated,  a  representative  sample  should 
have  included  this  region  as  well. 

The  north  central  area  of  the  province,  shown  as  region 
F,  is  another  mixed  farming  area.  Edmonton  and  Vermilion 
are  located  practically  on  the  border  between  region  E, 
and  region  F,  hence  they  may  reflect  the  weather  of  one 
area  as  much  as  the  other. 

The  Peace  River  district  is  relatively  new  and 
relatively  segregated  from  the  rest  of  the  province.  This 
area,  designated  as  region  G,  was  represented  by  both 
Beaverlodge  and  Fairviev. 

The  weather  stations  chosen,  as  well  as  being 
geographically  distributed  throughout  the  province,  appeared 
to  represent  the  major  farming  areas  of  the  province. 

Hence,  even  though  the  sample  was  chosen  out  of  necessity, 
the  province  was  as  '-ell  represented  as  may  be  expected 
from  a  small  sample. 


■ 
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3.6  Preparing  the  Data  for  Processing 


Having  decided  upon  the  data  required  and  the  stations 
from  which  data  .vas  to  be  used,  it  vas  possible  to  order 
cards.  A  summary  of  the  stations  used,  and  the  years  con¬ 
sidered  is  given  in  Table  9.  The  number  of  years  used 
depended  primarily  upon  the  number  of  years  for  vhich  cards 
were  available.  The  cards  for  Edmonton  and  Calgary  mere 
available  on  campus  for  1921  onward  and  this  vas  .the  period 
used  for  these  t o  stations.  The  meteorological  Branch 
suggested  cards  be  ordered  in  lots  of  five  years,  either 
in  2  -  6  or  7  -  l  combinations,  and  as  far  as  possible 
data  vas  obtained  back  to  1922. 


Table  9-  _ Cards  Put 


on  Tape 


Sta  t  ion 


Sta  t  ion  No  . 


Years 


Edmonton 


Vermilion 

Calgary 


3012195 

1921  - 

1937 

3012200 

1938  - 

1961 

3012205 

1963  - 

1965 

3012203 

1962 

3016800 

1947  - 

1966 

3031080 

1939  - 

1966 

3031085 

1921  - 

1931 

3031100 

1932  - 

1938 

3033860 

1938  - 

1966 

3033890 

1922  - 

1937 

Lethbridge 


-31- 

Station 

Station  No. 

Years 

Medicine  Hat 

5034480 

1922  -  1966 

Beaverlodge 

3070560 

1922  -1966 

Fair view 

3072520 

1932  -  1966 

except  1956 

In  total  there  vere  102,000  cards  to  be  handled.  In 
order  to  facilitate  computer  handling,  these  cards  had  to  be 
stacked  on  tape.  One  2400  foot  tape  was  required.  The 
stacking  process  required  about  tvo  hour’s  of  computing  time 
once  the  cards  vere  ready  for  stacking.  The  first  time 
stacking  was  attempted,  several  cards  with  invalid  characters 
were  detected  by  the  card  reader,  and  the  job  had  to  be 
stopped.  The  invalid  characters  were  found  to  be  multiple 
punches,  and  occurred  in  either  column  79,  or  column  80 . 

In  order  to  identify  the  invalid  cards,  all  the  cards  were 
put  through  a  card  sorter  set  to  reject  anything  by  an  X 
in  column  79*  In  some  cases  as  many  as  four  numbers  were 
punched  in  one  column.  Several  legal  characters  eg.  P,N, 
D,Q,$,&,  and  =  were  discovered,  in  addition  to  at  least  as 
many  different  illegal  characters.  This  might  indicate 
special  marking  since  they  vere  found  at  all  stations,  and 
only  between  1957  and  1962.  Altogether,  there  were  600  cards 
with  illegal  characters.  Since  the  illegal  characters 
occurred  only  in  the  last  two  columns,  the  cards  were 
duplicated  omitting  the  last  two  columns,  so  that  the 
information  was  not  lost.  The  cards  vere  then  stacked  on  tape. 
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Once  the  cards  vere  stacked  on  tape,  they  vere  sorted 
first  according  to  station  identification,  and  then  accord 
to  date.  This  means  that  because  a  station  number  did  not 
change  in  the  middle  of  s  year,  the  days  vithin  each  year 


ing 


are  in  chronological  order.  Hovever,  since  the  station 
numbers  did  not  necessarily  get  larger  when  they  vere 
changed,  the  years  are  not  necessarily  in  chronological 
order . 

The  tape  vas  then  checked  for  completeness  of  data. 
The  Computing  Science  Department  suggested  first  writing 
a  program  in  Cobol,  to  check  for  unexpected  data  in  a 
particular  column.  Since  the  writer  had  no  knowledge  of 
Cobol,  a  Fortran  program  vas  written  to  check  the  data 
column  by  column  for  the  columns  of  interest.  The  station 
number,  the  date,  minimum  and  maximum  temperature,  and 


total  precipitation  were  the  only  information  of  interest 
on  the  card  hence,  only  columns  1  -  19,  and  columns  45  -  48 
vere  checked.  Records  (or  cards)  with  anything  but  a  ;  in 
column  1,  a  0  in  column  2,  and  a  1,  5,  or  7  in  column  3 
were  written  out.  Columns  1  -  19  were  checked  for  anything 
but  numeric  data,  except  for  columns  14  and  17  which  could 
contain  a  blank,  a  zero,  or  a  minus  sign.  The  data  was 
also  checked  for  months  greater  than  12,  as  well  as  for  days 
greater  than  31 .  Columns  45  -  48  •••ere  checked  for  anything 


but  numeric  characters 


blanks,  or  the  letter  T. 


The  letter 


* 
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L  in  column  48  indicates  that  it  is  not  known  whether 
precipitation  occurred,  but  if  it  did  occur,  it  will  be 
included  in  the  next  day’s  total.  The  letter  C  in  column 
48  indicates  that  precipitation  did  occur  ori  that  day,  but  it 
will  be  included  in  the  next  day’s  total.  Although  the 
letters  L  and  C  were  expected,  it  was  of  interest  to  know 
where  they  occurred,  and  they  were  written  out  as  veil.  Table 
10  is  a  summary  of  days  with  C  or  L  in  column  48. 

A  second  program  was  used  to  locate  missing  days. 

This  program  had  limited  use.  First  the  data  had  to  be  in 

chronological  order  within  years.  This  presented  difficulties 

when  blanks  occurred  in  the  date  field  because,  for  example, 

6612  7  would  have  been  placed  before  661201  in  a  sort 

program.  This  caused  the  checking  program  to  lose  count  and 

write  out  every  day.  Alphabetic  data  could  not  be  read  in 

an  I  format,  and  an  error  condition  arose  when  alphabetic 

data  occurred  in  columns  1  -  13.  When  these  misfit  records 

were  located  it  was  possible  to  program  around  them.  This 

was,  however  a  hit  and  miss  procedure. 

Table  10  Days  with  Unknown  Precipitation 

Column  48 

Station  Dajte  C  L 


3012195  Aug  9,  1925  X 

Feb  22,  1939  X 

May  2,  19*10  X 

Jun  19,  19*10  X 


Sep  15,  1940 


X 
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Column  48 

Station 

Date 

C  L 

3031100 

May  19,  1931 

X 

Dec  4,  1924 

X 

May  28,  1930 

X 

Jul  18,  1930 

X 

Sep  15,  1931 

X 

3070560 

Aug  8,  1926 

X 

Aug  3,  1962 

X 

Sep  28,  1963 

X 

Feb  17,  1965 

X 

Sep  17,  1966 

X 

After  chec 

king  the  tape  with  both 

programs ,  mis sing 

data  were  found  to  occur  in  several  forms.  Some  days  were 
missing  completely,  some  days  had  missing  temperature, 
some  days  had  negative  precipitation,  and  others  had  been 
nispunched  making  them  unusable.  Records  with  negative 
precipitation  were  found  to  follow  records  containing  a  G 
or  an  L  in  column  48.  According  to  the  Meteorological 
Branch,  an  X  was  punched  in  column  45  to  indicate  a 
cumulative  precipitation.  However,  the  IBM/36 0  interpreted 
this  character  as  a  minus  sign,  which  resulted  in  a  neg¬ 
ative  precipitation.  Table  11  is  a  summary  of  missing 
data  indicating  what  data  were  missing,  and  what  was  done 


about  it. 

If  the  correction  is  not  specified,  then  the  data  were 


taken  from  the 


tat ion  listed  in  the  monthly  meteorological 
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27 

summaries.  There  were  88  missing  days,  46  missing 
temperatures,  16  days  vith  negative  precipitation,  and  10 
mispunched  cards  totalling  loO  days  vith  missing  data. 

This  amounted  to  0.086$  missing,  0.045$  vith  missing  temp¬ 
eratures,  0.016$  vith  negative  precipitation,  and  0.010$ 
mispunched,  or  0.157$  of  the  cards  vith  missing  data. 

A  comole te  taoe  was  created,  which  can  be  used  with 
programs  requiring,  minimum  and  maximum  temperature,  and 
total  precipitation.  Allowance  must  be  made  to  handle  a 
T,  C,  or  an  L  in  column  48  otherwise  the  tape  may  be  used 
without  difficulty. 
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4 .  TRACTABILI TY 


4.1  Estimating  the  Shear  Strength  of  Soils 

. . 29"; 

Nichols  (1932)  experimented  with  four  clay  soils,  a 
clay  loam,  and  a  silt  loam.  He  estimated  shear  strength 
experimentally  using  different  confining  pressures.  The 
apparatus  used  can  be  seen  in  Figure  3. 

Figure  3  Apparatus  for  Measuring  Shear  Strength 


centre  section  vas  moved  by 
lever  and  the  shear  force 
measured . 


Using  this  technique,  Nichols  found  the  maximum  shear 
force  at  or  near  the  plastic  limit  (near  maximum  cohesive 
force).  Maximum  shear  force  depends  upon  confining  pressure, 
as  well  as  upon  moisture  content.  A  graphical  example  of 
Ni chol’s  results  can  be  seen  in  Figure  4. 


Figure  4.  She ar  St re ng th  with 


Varying  Moisture  and  Confinement 


Pressure 


10 


80 
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Expressing  these  results  in  equation  form  shear 
strength  is  given  as: 

Fs  =  (Pu  -  M) (0 . 06Pn  +  P  +  1.6) 

Fn  • 


-••here  Fs  =  shear  strength 


Pu  «  upper  plastic  limit  (liquid  limit) 

M  ~  moisture  content  (fo) 

Pn  ~  plasticity  number  =  0.6C  -  12 
C  =  percent  clay  content 
P  =  pressure  in  psi. 

This  equation  is  not  accurate  for  soils  vith  a  low 
plasticity  number  (lov  clay  content).  Shear  resistance  in 
pure  sand  is  dependent  upon  particle  shape  (internal 
resistance),  and  confining  pressure,  and  shows  no  increase 
in  shear  strength  vith  variation  in  moisture  content.  The 
relationship  for  sand  may  be  expressed  as: 

Sf  “  P  tan0 

where  t anfl  -  coefficient  of  friction 
and  P  ~  confining  pressure. 


Silty  soils  on  the  other  hand,  are  quite  unpredictable, 
and  may  behave  any  here  between  a  sandy  soil  and  a  clay 
soil.  Because  of  the  difficulty  in  determining  shear 
resistance  in  non-plastic  soils  only  plastic  soils  have 
been  considered.  This  being  the  case  the  number  of  non- 
working  days  will  be  overestimated  for  non-plastic  soils. 
However,  this  is  perhaps  more  desirable  than  underestimating 


* 


- 


the  number  of  non-forking  flays. 

According  to  Nichols  the  shear  strength  chat  a  soil 
is  capable  of  developing  is  dependent  upon  confining 
pressure  and  moisture  content.  Figures  5,  6  and  7  illustrate 
the  confining  pressure  required  to  develop  a  particular 
shear  stress  for  a  range  of  moisture  contents  for  three 
particular  soils  in  the  province  of  Alberta.  A  description 
of  Alberta  soils  is  given  in  section  4,2.  The  shape  of 
the  curves  is  the  same  for  each  soil,  except  that  they 
shift  to  the  left. or  right  depending  upon  the  liquid  limit, 
and  they  shift  up  or  dovn  depending  upon  clay  content. 


Figure  5.  Shear  Strength  Curves  for  Malmo  SIC 


Figure 


6 


Shear 


Strength 


Cur 


ve; 


for 


Seven 


Persons  C 


Figure  7.  Shear  Strength  Curves  for  Coa'Ldale  CL-C 


Mo i s  t  ur e  Content 


The  confining  pressure  exerted  by  a  tractor  tire  is 
not  easily  determined,  but  if  the  resistance  of  the  tire 
is  assumed  to  be  neglibible,  the  pressure  can  not  exceed 
the  inflation  pressure  of  the  tire.  The  recommended  in¬ 
flation  pressure  for  larger  tractors  in  common  use  in  Alberta 
(eg.  the  John  Deere  5020,  or  the  Massey-Ferguson  1130)  is 
1 6  psi.  This  means  that  shear  strengths  requiring  a 
confining  pressure  in  excess  of  lopsi,  are  not  obtainable, 
making  the  part  of  the  curve  about  16  psi0,  of  no  use  in 
determining  the  shear  strength  available  to  pull  an 
implement.  If  at  a  given  moisture  content,  a  confining 
pressure  of  more  than  1 6  psi.  is  required  to  develop  the 
necessary  shear  strength  for  field  operation,  then  field 
operations  are  not  possible  at  that  moisture  content. 

4.2 _ S oil s  in  the  Province 


Because  the  shear  strength  of  a  soil  depends  upon  its 
clay  content,  as  veil  as  a  number  of  other  soil  properties, 
these  soil  properties  have  been  examined  for  the  province . 
Table  12  indicates  the  mechanical  analysis  of  surface  samples 
of  Alberta  soils,  as  determined  by  members  of  first  year 


soils  courses. 
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Table  12 

K  e  ch  a  n  i  c  a  1  A  na  1  y  s 

is  of  A 

lberta  Soils  (Ap  horizon) 

Soil 

liecha 

■  L 

Analysis 

Soil  Zone 

Soil  Series 

Class 

£s 

i 

go 

Brovn 

Cavendish 

SL 

77 

11 

12 

For e mo s  t 

L 

45 

35 

20 

Seven  persons 

C 

18 

32 

50 

Dark 

Carmangay 

LS 

81 

10 

9 

Brovn 

Gran  vim 

I, 

42 

37 

21 

Coaldale 

CL-C 

28 

30 

42 

Thin 

Black 

Irma 

LS-SL 

77 

12 

11 

Elnora 

L 

37 

38 

25 

Three  Hills 

C 

16 

36 

58 

Black 

peace  Hills 

LS 

82 

10 

8 

Angus  Ridge 

L 

42 

33 

25 

Mai  mo 

SiC 

13 

35 

50 

Grey 

Culp 

LS 

80 

8 

12 

Wooded 

Breton 

L 

42 

36 

22 

Mayv/ood  CL-SiCL 

30 

42 

28 

A  number  of  other  soil 

char a c 

ter is tics 

vere 

of  interest 

and  are  included  in  Tables 

13  -  17 

.  These 

soil 

properties 

-  ere  also 

determined  by  the 

members  of  the 

first 

year  soils 

courses . 


I.  in 
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Table  13 

L 1 qu id  L i m it  of 

Alberta  Soils  (A 

p  horizon) 

LS-SL 

L 

SiCl-C 

Brown 

19.1 

28.3 

45.5 

Dark  Brown 

19.1 

26.2 

42.7 

Thin  Black 

S3. 3 

39.3 

53.6 

Black 

19.8 

41.7 

59.7 

Grey  Wooded 

13.6 

25.9 

36.1 

Table  l4 

Field  Capacity 

of  Alberta 

Soils 

(Ap  horizon) 

LS-SL 

L 

SiCl-C 

Brown 

12.0 

21.8 

35.7 

Dark  Brown 

12.4 

23.9 

35.3 

Thin  Black 

1.6.3 

32.3 

46.5 

Black 

11.4 

33.0 

45.0 

Grey  Wooded 

9.6 

22.5 

29.3 

Table  15. 

Organic  Matter 

of  Alberta 

Soils 

(Ap  horizon) 

LS-SL 

L 

SiCl-C 

Brown 

2.6 

2.6 

4.9 

Dark  Brown 

1.9  . 

3.9 

3.2 

Thin  Black 

4.0 

8.3 

7.0 

Black 

2.6 

10.9 

12.8 

Grey  Wooded 

1.7 

2.7 

3.7 

T< 


Bulk  Density 

(A.S.G.)  of  A 

inert a 

Soils  (A] 

LS-SL 

L 

SiCl-C 

Brown 

1.4 

1.3 

1.1 

Dark  Brown 

1.4 

1.3 

1.2 

Thin  Black 

1.3 

1.2 

1.2 

Black 

1.4 

1.0 

1.1 

Grew  Wooded 

1.3 

1.2 

1.1 

Permanent  Wilting  Point  of 

Albert 

a  Soils  j 

LS-SL 

L 

SiCl-C 

Brown  . 

O'! 

• 

CD 

10.3 

22.3 

Dark  Brown 

5.2 

9.6 

18. 1 

Thin  Black 

7.2 

15  n 

23.5 

Black 

5-9 

l.e.5 

23.5 

Grey  Wooded 

3.7 

7.6 

11.1 

The  program  developed  by  Baler  and  Robertson  considered 
a  fixed  depth  of  moisture  rather  than  a  fixed  depth  of  soil. 


It  was  also  of  interest  to  knov  the  depth  of  soil  equivalent 
to  one  inch  of  available  water.  The  following  equation  may 
be  used  to  determine  soil  depth. 

Depth  of  water  =  <3 bn 

rou 

where  d  -  depth  of  the  soil  horizon 
b  ~  bulk  density  (A.S.G.) 
n  -  percent  available  vater 

Available  vater  may  be  found  by  subtracting  the  moisture 
content  at  the  permanent  wilting  point  from  the  moisture 


-4 


O 

O* 


conte 

nt  at  field  capacity. 

Table 

1.3  is  a  summary  of 

death 

of  soil  per  inch  of  v 

a  ter  for  Alberta 

Soils. 

Table 

18 .  Depth  of  Soil  Per  Inch 

of  Avails 

ble  Water 

LS-SL 

L 

SiCl-C 

Brown 

12” 

7" 

7" 

Dark  Brown 

10" 

r  M 

O 

5" 

Thin  Black 

8" 

3" 

4" 

Black 

15" 

6" 

4" 

Grey  Wooded 

12" 

6" 

r-  !! 

O 

4.5 

Determination  of  Soil 

Tr  .  c  t 

bill ty 

In  order  to  determine  the  moisture  content  above  which 

field  work  was  impossible,  a  John  Deere  5020  '-as  used  as  an 

example.  This  tractor  ivS  normally  equipped  vith  24.5:52 

rear  tires,  and  11.00:16  front  tires,  vith  outside  diameters 

28 

of  71  in.  and  58  in.  respectively.  The  Nebraska  Test  , 
done  vith  18.4:58  rear  tires,  indicates  the  JD  5020  has  a 

wheel  base  of  102  -  106  in.,  and  a  drawbar  height  of  20 

18 

in.  Traction  Tests  indicate  that  a  24.5:52  tire,  and  an 

18.4:58  tire  have  approximately  the  same  loaded  radius, 

vhich  means  the  drawbar  height  will  be  the  same  vith  either 

tire.  The  weight  distribution  vith  24.5:52  tires  (single) 

vith  liquid  ballast,  and  cast  iron  ballast  is  15,900  lbs. 

on  the  rear,  and  4800  lbs.  on  the  front. 

18 

Traction  Tests  indicate  that  the  maximum  drawbar  pull 


. 


' 


-49- 


is  42.5^  of  the  weight  on  the  rear  wheels,  vhich  in  this 

case  amounts  to  5900  lbs.  In  addition  to  drawbar  pull, 

consideration  must  be  given  to  the  effect  of  rolling 

resistance  on  the  required  shear  strength,  ho  work  was 

found  soecifically  relating  rolling  resistance  of  a  tire, 

25  7 

to  soil  moisture  content.  Janosi  and  Barger  et  al. 

indicate  that  rolling  resistance  is  a  function  of  tire  diameter 

Using  the  relationship  given  by  Barger  et  al . ,  the  coefficients 

of  rolling  resistance  (C)  for  tilled  loam  are, 

-0.8 

G  -  4.2(71)  -  0.088  (rear) 

-0.7 

G  =  4.2(53)  =  0.145  (front) 

Rolling  resistance  may  be  determined  by  multiplying 
these  coefficients  by  the  respective  weights  on  the  rear 
and  rront  wheels,  but  in  doing  this  •■•eight  transfer  must 
be  considered.  Weight  transfer  is  determined  by  taking 
moments  about  a  point  on  the  ground  directly  below  the  rear 
axle,  in  vhich  case  weight  transfer  is  equal  to  the  pull, 
times  the  drawbar  height,  divided  by  the  wheel  base,  or 
1150  lbs.  The  tractor  weight  is • now  distributed  15,000  lbs. 
on  the  rear,  and  5700  lbs.  on  the  front.  The  rolling  re¬ 
sistance  (R)  under  these  conditions  is, 

R  =  0.088  (15,000)  +  0.145  (5700)  »  i860  lbs. 

The  coefficients  of  rolling  resistance  (C),  and  rolling 

resistance  (R)  on  Fall  Rye  seeding  are, 

-0.7 

C  =  1.6  (71)  =  0.042  (rear) 


. 
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-0.7 

C  -  1.6  (38)  =  O.Oo^  (front) 

R  =  0.042  (15,000)  +  0.0c3  (5700)  *  860  Lbs. 

If  the  rolling  resistance  is  i860  lbs.,  and  the 

implement  draft  is  5900  lbs.,  the  total  draft  is  7760  lbs. 
Hcvever,  if  the  rolling  resistance  is  only  8c0  lbs.  che 
co cal  draft  is  6760  lbs.  Traction  Tests  indicate  that 

the  contact  area  of  the  rear  tires  vith  the  soil  is  384 
0**] »  in  • ,  v/ntcii  noans  tnut  on  tilled  iomn  the  recuired  she 


al 


strength  is  7760/884  =  8.8  psi.,  and  on  Pall  Rye  seeding 

6760/884  =  7.7  psi. 

34 

Reece  indicates  that  the  maximum  shear  value  chat  a 
soil  will  develop  is  obtained  only  at  100/  slip,  and  that 
at  a  normal  slip  of  15/,  only  80/  of  the  maximum  shear  is 
obtained.  This  means  that  the  maximum  soil  shear  strength 
required  under  the  above  conditions  is  11  psi.  and  9.6 
psi.  for  tilled  loam  and  fall  Rye  seeding  respectively. 

Assuming  a  confining  pressure  of  16  psi.,  a  shear 
strength  of  11  psi.  is  developed  at  49.4/  moisture  concent 
in  a  Maimo  silty  clay  -  clay  soil  (Figure  5),  *-’here  as  9.6 
psi.  shear  strength  is  developed  at  50. 7 Z  moisture.  Both 
these  moisture  contents  are  slightly  above  field  capacity 
(45/).  In  a  Seven  Persons  clay  soil  (Figure  6),  however, 

11  psi.  shear  strength  is  developed  at  35/  moisture  content, 
and  9.6  psi.  shear  strength  is  developed  at  36.3 Z  moisture. 
In  this  case  field  capacity  (35-7/)  is  between  these  t”o 
moisture  concents.  A  similar  situation  is  found  in  a 


.  . 
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Coaldale  clay  loam  -  clay  (Figure  7),  ’’here  a  shear 

strength  of  11  psi.  is  developed  at  3^  -S$  moisture,  and 

9.6  psi.  shear  strength  is  developed,  at  35-9/5  moisture. 

As  in  the  Seven  Persons  clay,  field  capacity  (35-3/5)  falls 

between  these  two  moisture  contents.  In  each  case  a  small 

moisture  change  makes  an  appreciable  change  in  shear  strength 

In  two  of  the  soils  described  above,  the  critical 

moisture  content  appears  to  be  in  the  neighborhood  of  field 

12 

capacity.  Browning  suggests  tha t  optimum  tillage  con¬ 


ditions  occur  at  a  moisture 


content  slightly  below 


field 


capacity.  Hence  for  these  tvo  soils  as  veil  as  for  non¬ 
plastic  soils,  it  seems  reasonable  to  choose  a  moisture 
content  slightly  below  field  capacity,  as  the  criterion  for 
classifying  a  day  as  a  work  day  or  a  non- work  day. 


In  the  third  soil  described  above  the 


cr  it  i ca  1  mo  i s  ture 


content  appears  to  be  above  field  capacity.  Buckman  and 

13  3 

Brady,  and  Baver  indicate  that  the  structure  of  fine 
textured  soils  may  be  damaged  if  tillage  operations  are 
performed  at  high  moisture  contents.  Hence,  good  farming 
practice  ’  ould  prohibit  tillage  at  moisture  contents  above 


field  capacity, 


and  the  choice  of  a 


critical  moisture  content 


slightly  below  field  capacity  will  not  give 


criterion 


that  is  too  low. 

fillability  of  fine ' textured  soils  may  be  determined  by 
the  optimum  moisture  for  good  soil  structure,  rather  than 
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vhether  or  not  it  is  possible  to  pull  s  tillage  implement 


i  u 


the  field. 


Hovever,  in  an  attempt  to  limit  the  number  of 


variables,  the  criterion  for  classifying  a  day  as  a  vcrk 
day  or  a  non-vork  day  :-as  based  as  far  as  possible  upon 


the  possibility  of  pulling  an  implement  in  the  field. 

In  a  clay  -  clay  loam,  one  inch  of  available  soil 

moisture  corresponds  to  approximately  five  inches  of  soil, 

or  the  too  three  zones  of  the  Soil  Moisture  Budget.  Severs 
‘  40 


and  Severs 

ind ica ce 

that  pressure  under  a 

found 

at  ion 

decreases  * 

ith  depth. 

By  a na 1 ogy  pressure 

under 

a  uire 

should  deer 

ease  vith 

depth,  making  soil  at 

lover 

depths 

important  in  determining  soil  tractability .  This  being 


the  case,  and 
Budget  did  not 
except  during 


the  fact  that  the  Versatile  Soil  Moisture 
indicate  moisture  removal  from  Zone  4, 
the  period  from  jointing  to  the  hard  dough 


stage  (June  17  -  September  3),  only  the  top  three  zones 
vere  used  to  classify  a  day  as  a  vcrk  day  or  a  non-vork  day. 
In  coarser  textured  soils  only  the  top  tvo  zones  needed  to 
be  considered,  since  one  Inch  of  Vater  vas  held  in  a  greater 


depth  of  soil . 


If  only  the  top  tvo  zones  are  considered,  it  is  not 
necessary  to  have  an  exact  estimate  of  the  moisture  content 
above  vhich  tillage  operations  are  not  possible.  The 
reason  for  this  is  that  small  moisture  changes  make  relatively 


large  differences  in  the  shear  strength  of  the  soil,  and  in 
the  top  tvo  zones  moisture  changes  are  quite  rapid.  Hence 


. 

the  time  required  for  the  soil  to  dry  sufficiently  is  quite 
short.  In  this  case  the  number  of  non- work  days  viLl  be 
less.  If  the  third  zone  is  considered,  then  exact  critical 
moisture  estimates  are  mere  important,  since  the  moisture 
change  is  much  s lover  and  a  1%  moisture  change  might  take 
one  or  to  o  days  . 

Since  a.  season  includes  the  months  April  through 
October  some  consideration  should  be  given  to  snowfall . 

The  Soil  Moisture  Budget  treats  snov  as  an  equivalent  amount 
of  rain.  However,  twelve  inches  of  snov  would  probably 
create  mere  non-work  days  than  a  one  inch  rainfall.  R.W, 
Longley,  Professor  of  Meteorology  suggested  that  if  there 
was  snow  on  the  ground  in  the  Spring  of  the  year,  the 
temperature  would  not  likely  rise  above  4-5°F.  In  the  Fall 
because  of  heat  in  the  ground,  the  snow  may  have  been  gone 
at  a  temperature  of  40°F.  Hence,  it  was  considered  that 
there  was  snov  on  the  ground  in  the  Spring  when  temperatures 
did  not  exceed  45°F.  during  the  day,  and  in  the  Fall  when 
temperatures  did  net  exceed  4o°F.. 


. 

' 

5  RESULTS 


5.1  Classifying  the  Data  as  Work  Days  and  Non-Work  Days 

The  moisture  content  above  which  it  was  impossible  to 
perform  tillage  operations  appeared  to  be  at  or  near  field 
capacity.  Since  field  capacity  *  as  the  upper  limit  of  the 
Soil  Moisture  Budget  a  moisture  content  of  99*5$  of 
capacity  seemed  to  be  a  reasonable  choice  of  criteria  for 
classifying  a  day  as  a  work  day  or  a  non-work  day.  A  trial 
run  was  made  using  criteria  of  99*5$  of  field  capacity  in 
the  top  three  zones,  with  no  consideration  for  snowfall. 

Based  on  the  results  of  check  data  (described  in  section 
5.9),  the  number  of  non-vork  days  appeared  to  be  under¬ 
estimated.  It  was  also  apparent  from  this  trial  run,  that 
some  consideration  of  snowfall  would  improve  the  correlation 
in  early  Spring  and  late  Fall. 

Since  this  trial  run  seemed  to  underestimate  the  number 
of  non-work  days,  a  second  set  of  conditions  was  used  to 
classify  the  data.  In  the  second  run  a  day  was  classed  as 
bad  if  the  moisture  content  in  any  one  of  the  top  three 
zones  was  above  95$  of  field  capacity,  or  if  there  was  snow 
on  the  ground.  The  inclusion  of  the  top  three  zones  was 
only  justified  in  medium  or  heavy  soils.  In  sandy  soils 
only  the  top  two  zones  needed  to  be  considered.  Hence  an 
additional  run  was  made  to  determine  probability  values  for 
sandy  soils.  In  this  run  the  criteria  were  the  same  except 
that  only  two  zones  were  considered. 


. 
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Arising  out  of  these  runs  two  classifications,  one  for 
heavy  and  medium  soil,  and  one  for  light  soil,  vers  developed 
and  are  shown  in  Appendix  8.1.  Both  of  these  runs  shoved 
much  closer  correlation  with  the  check  data,  than  did  the 
first  run. 


In  classifying  the  data,  a  work  day  was  designated 
as  an  X  and  a  non-work  day  was  designated  as  a  (.).  In 
order  to  make  the  X!s  and  dots  more  useful  summaries  were 


made  of  the  frequency  of  runs  of  n  bad  days,  frequency  of 
runs  of  n  good  days,  and  the  number  of  years  with  n  bad 
days  in  a  month.  These  summaries  have  been  included  in 


Tables  19  -  24. 
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Table  21  Number  of  Years  with  n  Bad  Days  Per  Month  (Heavy  to  Medium  Soil) 
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5_ ._2_  Persia  tence  of  Goo d  and  Bad  Days 

Using  the  information  pertaining  to  the  frequency 
of  runs  of  n  bad  days,  and  of  n  good  days,  an  estimate  of  the 
average  number  of  bad  and  good  days  "as  made  for  each  month. 
Knowing  the  average  number  of  bad  and  good  days,  it  "as  then 
possible  to  determine  the  probability  of  a  day  being  either 
good  or  bad  in  any  given  month.  The  probability  of  a  day 
being  good  or  bad  varies  depending  upon  the  previous  day's 
weather . 

In  order  to  examine  this  aspect  more  closely  the  month 

of  May  at  Edmonton  may  be  considered.  Over  a  period  of  45 

years,  1920  -  1965  there  vas  an  estimated  total  of  363  bad 

days  in  medium  to  heavy  soil.  This  means  that  the  general 

probability  of  a  bad  day  is  0.262.  If  the  probability  of 

a  day  being  bad  is  designated  as  p,  and  the  probability  of  a 

day  being  good  is  designated  as  q,  then  the  probability  that 

a  day  begins  a  run  of  bad  days  is  pq.  If  the  chance  of  the 

following  day  being  bad  is  also  p,  then  the  probability  of  a 

2 

run  of  at  least  t* o  days  is  qp  ,  of  at  least  three  days 

etc.  providing  there  is  no  "persistence".  That  .is, 
if  the  probability  of  a  day  being  bad  is  not  greater  than  p 
when  the  previous  day  vas  also  bad.  The  total  number  of 
days  is  L395,  hence  the  expected  number  of  runs  pqN  is 
(1395)  (0.262)  (0.738)  -  270.  The  expected  number  of  runs 
of  at  least  two  bad  days  is  (270)  (0.262)  =  70.6,  of  at  least 


•1  T  ?  .  ,.*5i  ‘  .!J'  •  (  "  V-  '  t 
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three  bad  days  (70.6)  (0.262)  =  18.5  and  so  on.  These  values 
are  given  in  the  first  line  of  Table  25.  The  second  line 
indicates  the  number  of  runs  of  exactly  n  bad  days  vhich  is 
found  by  subtracting  successive  values  in  line  one.  The 
third  line  indicates  the  observed  number  of  runs  of  n  bad 
days . 


Ta b 1 e2 5  .  Fr e q uency  of  Runs  of  Bad  Da y s  a t dmonton  in  May , 


19^2  -  1965 

- — — - 

*  ■'  J  i 

n 

1 

2 

5 

4 

5 

6 

7 

8 

Q 

y 

10 

Runs  of  n  or 
more  days 

270. 

70.6 

18.5 

4.8 

1.3 

0.3 

0.1 

0.0 

0.0 

0.0 

Runs  of  n  days 

199.4 

52.1 

13.7 

3.5 

1.0 

0.2 

0.1 

0.0 

0.0 

0.0 

Observed  runs 

43. o 

16.0 

12.0  11 .0 

9.0 

4.0 

7.0 

7.0 

3.0 

!  .0 

It  vould  appear  that  there  are  too  many  one  day  runs, 


and  too  few  longer  runs  suggesting  possible  persistence. 

11 

Brooks  and  Carruthers  examined  persistence  of  rainy  days 


at  Kev,  and  suggest  that  there  are  three  possibilities 
concerning  persistence  probabilities. 

i  the  probability  of  rain  on  a  day  following  a  rainy 
day  is  constant  (i.e.  pk  =  P]_  for  all  values  of  k) 
ii  the  probability  of  rain  after  two  rainy  days  is 


constant  (i.e.  p1<;  =  p^  for  all  values  of  k  greater 
than  l) 


iii  p^  continually  increases  as  k  increases. 


Consider  the  first  possibility  first. 


The  f r e q ue n cy  o f 


2 

> 


5, 


observed  runs  of  l, 


,  n  or  more  days  is 
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113  (1,  Pp,  Pp^,  . ,  P]_n)  the  sum  of  T-hich  is  113/ 

(1  -  p^).  This  sum  may  be  equated  to  the  observed  number 
of  bad  days  363  yielding  -  0.688.  The  calculated  sei’ies 
of  runs  of  n  days  can  nov  be  seen  in  Table  29* 


Table 

26 

Freque; 

ncy  < 

of 

Runs 

of 

n 

Bad 

D 

ays  a 

t  Edmonton 

n 

1 

2 

3 

4 

5 

6 

7 

8 

o 

> 

10 

Runs 

of  n 

or 

113 

.0 

78 

.0 

53. 

.8 

37* 

1 

25* 

,  6 

17 

.6 

12 

.2  8. 

4 

5. 

.8 

4.0 

more 

days 

Runs 

of  n 

days 

35 

.0 

24 

.2 

16, 

.7 

11 . 

5 

8, 

,0 

5. 

.4 

3 

.6  2. 

£ 

V J 

1 , 

.8 

1.2 

Observed  i 

.urns 

43 

.0 

16 

.0 

12, 

11. 

0__ 

9j 

4 

i  i 

,_0_ 

7 

.0  7* 

0_ 

Jj 

.0 

1.0 

The  calculated  frequencies  are  now  much  closer  to  the 
observed  frequencies.  In  fact,  a  Chi-Square  test  with  five 

p 

degrees  of  freedom  (x  =  7.17)  indicates  that  a  verse  fit 
will  occur  due  to  chance  alone  22  times  out  of  100,  A  fit 
as  good  as  this  makes  it  unnecessary  to  check  other  possible 
probability  theories.  Hovever  as  a  matter  of  Interest  vas 
calculated  and  found  to  be  0.395*  A  Chi  Square  test 
(x<  ~  2.74)  t hen  allowing  for  the  one  less  degree  of  freedom 
does  improve  the  fit. 

The  probability  of  a  bad  day,  given  that  the  previous 
day  was  good,  and  the  probability  of  a  bad  day,  given  that 
the  previous  day  was  bad,  was  calculated  for  each  month  at 
each  station.  This  was  also  done  for  good  days.  These 


results  are  summarized  in  Tables  27  -  30.  The  probabilities 
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in  parentheses  indicate  the  probability  of  a  bad  day,  given 
that  the  previous  day  vas  bad,  and  the  probability  of  a  good 
day,  given  that  the  previous  day  vas  good,  respectively. 

Except  for  the  months  of  April,  September,  and  October 
a  Chi  Square  test  indicated,  a  good  fit  when  the  probability 
of  a  bad  day  vas  assumed  to  be  p,  .  In  total  there  -.-ere 
fifteen  instances  for  which  the  Chi  Square  test  indicated  a 
poor  fit  at  the  5 %  level.  At  the  1  fo  level,  the  Chi  Square 
test  indicated  a  good  fit  for  eight  of  the  fifteen.  It 
seems  reasonable  then,  to  assume  that  the  probability  of  a 
bad  day  is  constant  after  having  one  bad  day.  This  vas  also 


the  case  with  good  days. 

If  one  wishes  to  simulate 


month’s  weather,  he  must 


make  some  allowance  for  persistence  of  weather.  In  the  month 
of  May  at  Edmonton  .for  example,  the  first  day  may  be 
classified  as  a  work  day  or  a  non-work  day  on  the  basis  of  the 
general  probability  0.26  for  a  bad  day,  and  0.7^  for  a  good 
day.  However,  having  classified  the  day  one  way  or  the 
ether,  the  general  probabilities  0.2c  and  0.7^  a r o  no  longer 
valid.  Instead,  if  for  example  the  first  day  vas  bad,  the 
probability  of  May  2  being  bad  is  Pj  cr  0.69.  The  probability 
of  May  2  being  good  is  1.0  -  O.69  -  0,31.  If  may  2  is  good 


then  the  probability  of  May  3  being  good  is  0.87.  The 
probability  of  May  3  being  bad  is  1.0  -  0.87  -  0.13. 

The  probabilities  p  and  can  be  used  to  estimate 


the 


. 
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theoretical  length  of  run  of  bad  days  assuming  no  per¬ 
sistence  a nd  ass urn i ng  persiste n c e  respectively.  A s s urn i ng 
no  persistence,  the  length  of  run  vould  be  l/(l-p)  - 
l/(l.O  -  0.2c)  -  1.35  days.  Assuming  persistence  does  occur, 
the  theoretical  length  of  run  of  bad  days  vould  be  i/U-Pi) 

=  1/(1. 0  -  O.Cg)  -  3-22  days. 

Persistence  of  vork  and  non- work  days  v'as  found  to 
occur  regardless  of  the  method  used  for  classifying  days 
as  vork  or  non-work  days.  Hence  the  need  for  making  allovance 
for  persistence  of  veather  events  is  clearly  pointed  out. 


5  .J5 _ A  t_t  snip  ted  Checks  o  n  Program  A  ccura  cy 

Since  the  relationship  between  soil  moisture  content  and 
soil  trac lability  has  not  been  veil  established,  it  vas 
desirable  to  have  some  assurance  that  the  estimates  of  good 
and  bad  days  were  as  accurate  as  the  Soil  Moisture  Budget. 

Tvo  methods  were  used  to  check  the  results.  The  first  type 
of  check  involved  a  comparison  between  observed  and  est¬ 
imated  working  days  when  precipitation  was  known.  In  the 
second  check  the  working  days  were  known,  but  precipitation 
was  not,  hence  a  direct  comparison  could  not  be  made. 

The  first  type  of  check  vas  limited  to  some  extent, 
in  that  only  one  year’s  data  were  available  for  two  stations. 
One  station  1  as  located  at  Rutledge  Ranching  Co.,  Monitor, 
Alberta,  and  the  second  station  was  located  at  Liivam  Farms 
Ltd.,  Penhold,  Alberta.  At  Monitor  a  record  vas  kept  of 
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precipita tion,  as 
on  a  daily  basis 
precipitation  vas 
the  soil  was  dry 


veil  as  any  field  operation  in  progress 
(Appendix  8.2.1).  At  Penhold  a  record  of 
kept,  and  an  estimate  vs s  made  as  to  whether 
enough  to  permit  field  operations 


(Appendix  8.2.2).  The  soil  at  Monitor  vas  a 
•fhere  as  the  soil  at  Penhold  vas  a  silt  loam. 


sandy  loam. 
Using  the 


first 


method  of  estimating  bad  days,  that  is  using  moisture  criteria 
of  99.5/£  of  capacity  in  the  top  three  moisture  zones,  the 
number  of  bad  days  were  underestimated  by  about  one  day 
for  one  inch  rainfalls,  and  by  two  to  three  days  for  two 
i  nch  r  a  i  nf  a  1 1  s  .  Us  ing  the  s  e  co  rid  me  the  d  of  est  ima  ting  bad 
days,  that  is  using  moisture  criteria  of  95/3  of  capacity  in 


the  top  three  zones,  at  Monitor  a  one  inch  rainfall  caused 
the  number  of  bad  days  to  be  overestimated  by  one  to  tvo 


days.  At  Fenhold  the  number  of  bad  days  were  underestimated 
about  the  same  number  of  times  as  they  were  overestimated. 

The  results  of  this  check  can  be  explained  by  examining 
the  soils.  At  Monitor  the  soil  vas  a  sandy  loam,  hence  a 


criteria  designed  for  finer  textured  soils  vould  over¬ 
estimate  the  number  of  bad  days.  The  soil  at  Penhold  vas 


finer  textured  and  estimates  of  bad  days  were  closer. 


Although  the  results  of  one  yea  ts  data  are  obviously  inconclus 
ive,  the  accuracy  of  the  program  appears  to  be  fairly  good. 

The  second  check  vas  made  possible  as  the  result  of  a 


coincidental  conversation  vith  Mr.  Don  Matthevs,  owner  of 
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the  Highland  Stock  Farm,  located  about 

Calgary  on  the  south  side  of  the  Elbow 

26 


ten  miles  vest  of 
River.  Mr .  Ma t the v s 


has 

and 


kept  a  diary  of  farm  operations  for  a  number  of  years, 
vas  kind  enough  to  make  the  diary  available  for  purposes 


of  checking  the  method  used  in  estimating  vork  days.  A 
fairly  complete  record  of  field  operations  vas  obtained  for 
the  years  1944  -  1966.  Although  the  amount  of  precipitation 
is  not  generally  recorded,  interruptions  in  field  vork  due  to 
rainfall  were  noted.  These  records  provided  an  alternative 
method  of  estimating  good,  and  bad  days. 

Although  the  farm  is  located  only  ten  miles  from 
the  city  limits,  it  is  located  twenty  miles  from  .'the  Calgary 
weather  station.  The  distance  involved  made  it  impossible 
to  assume  identical  precipitat ion  at  the  two  locations.  It 
is  possible  however,  that  over  a  period  of  twenty-three 
years  the  differences  in  precipitation  would  average  out. 

If  this  is  in  fact  true,  then  the  probability  of  a  bad  day 
at  one  station  should  be  approximately  the  same  as  the 
probability  of  a  bad  day  at  the  other  station. 

Tvo  methods  were  used  in  determining  the  probability 
of  a  bad  day  in  some  months.  Using  the  first  method,  if 
tillage  operations  were  in  progress,  or  if  it  appeared  that 
tillage  operations  could  have  been  carried  out  the  day  was 
classed  as  good.  If  field  operations  were  stopped  because 
of  rain,  or  if  it  appeared  that  it- vas  too  vet  to  vork  a 
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day  was 
when  no 
esti 


classed  as  bad.  However  in  some  instances 
tillage  operations  were  in  progress,  hence 
e  was  made  by  considering  only  the  days  whan 


it  rained 
a  second 
til lage 


operations  were  in  progress,  and  the  days  when  tillage 
operations  appeared  to  be  stopped  due  to  rain.  The 


probabilities  using  these  two  methods. 


as  well  as  the 


probabilities  determined  for  medium  to  heavy,  and  sandy  soil 
at  Calgary  are  given  in  Table  31* 


Table  31 

Probability 

of  a 

Bad  Da 

7  -  ch 

eck 

result 

s 

Apr 

May 

June 

July 

Aug 

Sept 

0  c  t 

Ma t  thews 

1  .79 

.26 

.21 

.10 

.12 

.17 

.22 

Ma  t thews 

O 

- 

.29 

.14 

.22 

- 

-- 

Calgary  ( 

'heavy)  .64 

.29 

.35 

.20 

.22 

•  36 

.51 

Calgary  ( 

[sandy)  .  4l 

.18 

.24 

.16 

.17 

.22 

.26 

6  CONCLUSION 


The 

weather,  v 


available 


object  of  this  study  was  to  examine  the  daily 
ith  a  view  to  establishing  the  number  of  days 
in  each  month  for  tillage  operations.  Criteria 


were  developed  to  classify  a  day  as  a  work  day  or  a  non- 
work  day.  Using  these  criteria  and  several  year’s  weather 
data,  the  probabilities  of  a  good  and  of  a  bad  day  were  d 
termined  on  a  monthly  basis  for  seven  Alberta  stations. 

One  obvious  conclusion  that  can  be  drawn  from  the 

results  is  that  persistence  of  weather  events,  a  ohenomenon 

11 

observed  by  Brooks  and  Carruthers  ,  creates  persistence 
of  work  days,  or  non-work  days  as  the  case  may  be.  This  means 
that  the  probability  of  a  day  being  either  a  work  day,  or  a 
non-work  day  can  not  be  established  objectively,  without  some 


knowledge  of  conditions  existing  on  the  previous  day.  That 
is,  the  probability  of  a  day  being  bad  when  the  previous  day 
was  good,  is  not  the  same  as  the  probability  of  a  bad  day 
when  the  previous  day  was  also  bad.  Similarity,  the  probability 
of  a  good  day  is  dependent  upon  the  conditions  of  the  previous 
day.  This,  of  course  lengthens  the  runs  of  good  and  bad 
days . 

The  probability  values  determined  using  criteria  of 
99.5/3  of  moisture  capacity  in  the  top  three  moisture  zones 
yielded  almost  the  same  values  as  those  determined  using 
criteria  of  95/3  of  moisture  capacity  in  the  top  two  zones 
except  in  the  months  in  which  snow  was  a  factor. 


This  may 


■ 


I  I 


_  P  q  _ 
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be  explained  in  tvo  vays .  In  the  case  of  a  large  rainfall, 
the  top  two  zones  are  dried  belov  95 $  of  capacity  by  the 
time  the  third  zone  is  dried  to  99*5$  of  capacity.  In  the 
case  of  smaller  rains  that  do  not  affect  the  third  zone,  it 
does  not  make  a  great  deal  of  difference  vhether  95 $  or 
99. 5^6  is  chosen  os  criteria,  since  the  time  required  to  reach 
either  moisture  content  is  quite  short. 

Since  tillage  operations  may  be  performed  in  coarse 
textured  soils  at  moisture  contents  near  field  capacity  vith 
out  harm  to  soil  structure,  the  probabilities  estimated  using 
either  method  should  satisfactorily  reflect  real  conditions. 
In  finer  textured  soils  tillage  at  high  moisture  contents  is 
harmful  to  soil  structure,  and  the  second  method  must  be 
used.  That  is,  tillage  operations  are  not  performed  until 
moisture  content  is  beloT'’  95$  of  capacity  in  the  top  three 
moisture  zones.  The  probability  of  a  bad  day  estimated  on 
this  basis  may  in  fact  still  be  conservative,  however  on  the 
basis  of  physical  ability  to  pull  a  tillage  implement,  the 
estimates  made  using  a  criterion  of  95$  in  the  bop  three 
zones,  appear  to  be  in  the  right  neighborhood.  Caution  must 
hovever,  be  exercised  vhen  using  these  probabilities  for  fine 
textured  soils,  like  Malmo  silty  clay  for  example. 

Although  the  probabilities  determined  here  may  not  be 


the  best  under  all  soil  and  topography  conditions,  they  do 
help  quantify  variables  of  great  significance  in  the 


. 
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scientific  business  of  farming  today.  Improvements  in  th- 
estimated  probabilities  for  tillage  operations  as  veil  as 
the  establishment  of  probability  values  for  harvesting 
operations  vill  no  doubt  provide  the  farmer  vith  a  better 
guide  to  correct  machinery  sizing. 
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5.  APPENDICES 


8.1  Listing  of  Work  and  Non  -  Work  Days 


8,1.1 _ Heavy  to  Medium  Soils 

In  making  this  estimate  of  good  and.  bad  days,  a  day 
vas  classified  as  bad  if  the  moisture  content  in  either  of 
the  top  three  soil  moisture  zones  vas  greater  than  95/j  of 
the  capacity  of  that  respective  zone,  provided  that  there 
vas  no  snov  on  the  ground.  It  vas  assumed  that  if  there 
vas  snov  on  the  ground  In  the  Spring  of  the  year,  the  air 
temperature  did  not  rise  above  45°F.  If  there  vas  snow  on 
the  ground  in  the  Fall,  it  vas  assumed  that  the  air  temperature 
did  not  rise  above  40°P,  Hence,  it  vas  considered  that  there 
vas  snov  on  the  ground  in  the  Spring. when  temperatures  did 
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8 . 1  M2 _ S  andy  Sol  Is 

In  making  estimates  of  good  and  bad  days  for  sandy 
soils,  a  day  vas  classified  as  bad  if  the  moisture  content 
in  either  of  the  top  t'  o  soil  moisture  zones  vas  greater 
than  95$  of  the  capacity  of  that  respective  zone,  provided 
that  there  vas  no  snow  on  the  ground.  Again  it  vas  assumed 
that  if  there  vas  snov  on  the  ground  in  the  Spring,  the 
temperature  did  not  rise  above  45°F.  If  there  vas  snov  on 
the  ground  in  the  Fall,  it  vas  assumed  that  the  air  tempera tur 
did  not  rise  above  4o°F.  Hence,  it  vas  considered  that 
there  vas  snov  on  the  ground  in  the  Spring  when  temperatures 
did  not  exceed  4p°F.  during  the  day,  and  in  the  Fall  vhen 
temperatures  did  not  exceed  4o°f. 


The  vea trier  stations  are: 
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8.8  Weather  Data  Used  .for  Checking  Purposes 


8.8.1  Work  Data  From  Rutledge  Ranching  Co. 


April,  1966 
Day 


Coera t ion 


Comments 


1 

8 

■z; 

4 

5 
b 

7 

8 

9 

10 
11 
18 
17 

14 

15 

16 

17 

18 
19 
80 
81 
88 

23 

24 

25 

26 

27 

28 

29 


Started  cultivating  vith  Noble  Blade 


Blading,  started  seeding  grass 


. 


May,  i960 


■S2~ 


D; 


Operation 


Comma  nt 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

1 7 

18 

19 

20 
21 
22 

23 

24 
28 

26 

27 

28 


Bladed  1100  acres ,  seeded  another 
1100  acres  since  starsing. 


No  rain  since  started 
field  vorlc 


29 

30 

31 


Thu nd e r  s ho v e r 
Rain  l/2  in. 


1/2  in. 


Ju 

Da; 

1 

2 

3 

4 

5 

6 

7 

8 

Q 

10 

11 

12 

13 

14 

15 

l6 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 
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1966 


Operation 


Comments 


Rain  1  in. 

Due  to  drought  conditions  seeding  of  cereal  crops 
had  not  started 

Dry  enough  to  plov  furrow  to  plant  trees 

Rain  2/10  in. 


Began  seeding  3:50  P.M 
soil,  other vise  okay 
Seeding 
Seeding 


soil  stuck  to  disc  in  heavy 


Seed ing 


No  field  work 
No  field  work 

Seeding  and  rodveeding  1:00 

Rod weed ing 

Seed ing  2 : 00  P.M . 


P.M. 


Cul t i va t ing  in  A . M . , 
Cultivating 
S e  e  d  i  ng ,  r  o  d  w  e  e  d  i  ng 
Seeding,  rodveeding 
Seeding,  rodveeding 
Seeding,  rodveeding 
S  e  e  d  i  ng ,  r  o  d r.  •  e  e  d  i  ng 
Seeding,  rodveeding 
Seeding,  r o  d  v e  e  d i ng 


seed ing 


Rain  l/lO  in*  after  supper 
Soil  damp  but  still  able 
to  work 
Rain  1  in. 


Rain  L/l0  in.  early  A.M., 
Rain  3/10  in.  6:00  P.M. 
Rain  1/10  in.  2:00  P.M. 


in  P.M. 


Rain  3/10  early  A.M. 

3/10  in  7:00  P.M. 

Good  weather 

Good  weather 

Good  weather 

Good  weather 

Good  weather 


X 
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July,  i960 


Day 

Operation 

v. 

Comments 

1 

No  field  vork 

2 

No  field  vork 

3 

Iio  field  vork 

Rain  8/10 

in. 

4 

No  field  vork 

5 

No  field  vork 

6 

No  field  vork 

7 

No  field  vork 

8 

No  field  vork 

Rain  2/10 

in. 

9 

No  field  vork 

10 

No  field  vork 

Rain  7/10 

in. 

ll 

Began  sva thing  hay  5:00  P.M. 

12 

Sva thing 

13 

No  field  vork 

Rain  2/10 

in  early  A.M. 

14 

No  field  vork 

15 

Sva  thing ,  cul t iva t ing 

16 

Sv a th ing ,  baling  5:00 

P.M. 

17 

Sv a  thing ,  bailing 

18 

Swathing,  baling 

19 

Baling 

Shower  in 

P.M.  stopped  vork 

20 

Baling  4:00  -  5:00  P. 

M.  Shower  in 

A.M.,  shower 

5:00  P.M. 

21 

Baling  10:00  -  1:00  A 

.M. 

22 

Baling  1:00  P.M. 

23 

Baling  8:00  A.M. 

Shower  5: 

00  P.M. 

24 

Baling  10:00  A.M.  **7:00  P.M.  (breakdown) 

7/10  in.  in  P.M. 

25 

Sva  thi. ag ,  bal  ing  8:30 

A.M. -11:00  A.M.  Rain  ' 

26 

No  field  vork 

27 

Baling  2:45  -  7:00  P. 

M. 

28 

Baling  7:30  A.M.  -  J: 

00  P.M. 

29 

Baling  7:30  A.M.  -  J: 

00  P.M. 

30 

No  field  vork 

Rain  5/10 

in.  early  A.M. 

31 

No  field  vork 

Rain  4/10 

in.  early  A.M. 

■ 


Operation 
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IQoo 


Comment! 


1 

No  field  vork 

2 

Baled  8:00  A.M.  ~  7:00  P.M. 

3 

Baled  8:30  A.M.  -  3:30  P.M. 

(breakdovn )  Ri in  .35 

4 

No  field  vork 

Rain  2/10  in. 

5. 

No  field  vork 

Rain  2  in. 

6 

No  field  Tork 

Rain  6/10  in. 

7 

No  field  vork 

8 

No  field  vork 

a 

Jr 

No  field  vorlc 

10 

Sva  th i ng ,  e ul  t  i va  t  i ng 

11 

Baling  5:00  -  6:15  P.M. 

Cloudy  all  day 

at  6:15  P.M. 

12 

Baling  10:30  A.M.  -  7:00  P.M 

• 

13 

Baling  8:00  A.M.  -  7:00  P.M. 

14 

Baling  8:00  A.M.  -  7:00  P.M. 

15 

Baling  8:00  A.M.  -  7:00  P.M. 

16 

Baling  8:00  A.M.  -  7:00  P.M. 

17 

Baling  8:00  A.M.  -  7:00  P.M, 

18 

Baling  1:30  -  7:00  P.M.  (tough)  Shovers  in  A.M 

19 

B'- led  9:00  A.M,  -  7‘?00  P.M. 

20 

Baled  8:30  A.M.  -  7:00  P.M. 

21 

Baled  8:30  A.M.  -  7:00  P.M. 

22 

Baled  8:30  A.M.  -  3:00  P.M. 

f breakdovn) 

23 

Baled  9:30  A.M,  -  7:00  P.M. 

(1  hr .  b r e a kdo v n ) 

24 

Baled  8:30  A.M.  -  7:00  P.M. 

25 

Baled  8:30  A.M.  -  7:00  P.M. 

26 

Broke  dovn  all  day 

27 

Finished  baling 

28 

29 

30 

31 

in. 

2:30 


6:30 

?.M. 


.  i'i 


ri  r 

-  - 


2 


Work 


Da 


f-  r* 


from  Li i van.  Farms  Ltd 


June,  1966 

Da y  Opera  bio n  Common t s 

1 

2 

3 

4 

5 

6 

7 

3 

9 

10 
11 
12 

13 

14 


15 

Field 

operations 

possible 

16 

Field 

operations 

possible 

1 7 

Field 

opera t ions 

possible 

18 

Field 

opera tions 

possible 

19 

Field 

opera  t ions 

possible 

20 

Field 

operations 

possible 

21 

Field 

ooera tions 

X 

possible 

22 

Field 

opera tions 

possible 

23 

Field 

operations 

possible  l/lO 

in. 

rain 

24 

Field 

opera  cions 

not  possible  2/10 

in . 

rain 

25 

Field 

opera tions 

no  t  possible  2/l 0 

in. 

rain 

26 

Field 

operations 

not  possible  1/10 

in. 

ra  in 

27 

Field 

operations 

possible  at  4:00  P.M. 

23 

Field 

opera t ions 

possible  at  11:00  A.M. 

1/10  in.  rain 

29 

Field 

operations 

possible 

30 

Field 

operations 

not  possible  3/10 

in. 

rain 

July 

,  I9c6 

Day 

Operation 

Comments 

1 

Field 

opera  tion 

not  possible 

1/10  in. 

ra  in 

2 

Field 

opera  t ion 

not  possible 

1  in.  rain 

3 

Field 

operation 

not  possible 

2.3  in. 

rain 

4 

Field 

operation 

not  possible 

1/10  in. 

ra  in 

5 

Field 

operation 

not  possible 

2/10  in. 

rain 

6 

Field 

operation 

not  possible 

7 

Field 

opera tion 

possible 

3 

Field 

opera tion 

possible 

9 

Field 

operation 

possible 

10 

Field 

operation 

possible 

11 

Field 

operation 

possible 

12 

Field 

operation 

possible 

13 

Field 

operation 

possible 

14 

Field 

operation 

possible 

15 

Field 

opera  tion 

possible 

16 

Field 

operation 

possible 

17 

Field 

opera tion 

possible 

18 

Field 

operation 

possible 

19 

Field 

opera  tion 

possible 

1/10  in. 

20 

Field 

operation 

possible 

ra  in 

21 

Field 

operation 

possible 

22 

Field 

operation 

possible 

23 

Field 

operation 

possible 

24 

Field 

operation 

possible 

6/10  in. 

25 

Field 

opera  tion 

not  possible 

rain 

26 

Field 

operation 

not  possible 

3/10  in. 

ra  in 

27 

Field 

opera  tion 

not  oossible 

28 

Field 

operation 

11:00  A.M. 

29 

Field 

operation 

possible 

30 

Field 

operation 

not  possible 

1.1  in. 

rain 

31 

Field 

operation 

not  possible 

August,  l$C6 


Day 

Operation 

Comments 

1 

Field 

operations 

possible 

2 

Field 

operations 

possible 

3 

Field 

operations 

possible 

0.1 

in.  rain 

4 

Field 

operations 

po  s s ib 1 e 

0.1 

in.  rain 

5 

Field 

operations 

not  possible 

1.4 

in.  rain 

6 

Field 

operations 

not  possible 

1.4 

in.  rain 

~7 

7 

Field 

operations 

not  possible 

0.2 

in.  rain 

8 

Field 

operations 

not  possible 

9 

Field 

operations 

possible 

10 

Field 

operations 

possible 

11 

Field 

operations 

possible 

12 

Field 

opera t ions 

possible 

13 

Field 

operations 

possible 

14 

Field 

operations 

not  possible 

0.3 

in.  rain 

13 

Field 

operations 

not  possible 

16 

Field 

opera  lions 

net  possible 

0.3 

in.  rain 

17 

Field 

operations 

not  possible 

18 

Field 

operations 

not  possible 

19 

Field 

operations 

net  possible 

0.2 

in.  rain 

20 

Field 

opera  tic ns 

not  possible 

0.1 

in.  rain 

21 

Field 

operations 

not  possible 

22 

Field 

operations 

possible 

23 

Field 

operations 

possible 

24 

Field 

operations 

possible 

25 

Field 

opera  t ions 

possible 

26 

Field 

opera t ions 

possible 

27 

Field 

operations 

po s s ib 1 e 

28 

Field 

operations 

not  possible 

1.0 

in.  rain 

29 

Field 

operations 

not  possible 

30 

Field 

operations 

not  possible 

0.2 

in.  rain 

31 

Field 

operations 

not  possible 

.. 
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Sept 

ember , 

196  c 

Day 

Opera t ion 

Comrne  nt 

3 

1 

Field 

opera  tions 

po s s lb 1 e 

2 

Field 

opera tic ns 

possible 

3 

Field 

operations 

possible 

4 

Field 

operations 

possible 

5 

Field 

operations 

possible 

6 

Field 

operations 

possible 

7 

Field 

operations 

not  possible 

0.2 

in. 

rain 

o 

Field 

opera t ions 

not  possible 

0.1 

in. 

rain 

9 

Field 

operations 

not  possible 

10 

Field 

operations 

possible 

11 

Field 

opera  t ions 

possible 

12 

Field 

operations 

possible 

13 

Field 

operations 

possible 

14 

Field 

operations 

possible 

15 

Field 

operations 

possible 

16 

Field 

operations 

possible 

17 

Field 

opera  t ions 

possible 

18 

Field 

opera  tions 

possible 

19 

Field 

operations 

possible 

20 

Field 

operatic ns 

possible 

21 

Field 

opera tions 

possible 

22 

Field 

00 era tions 

possible 

23 

Field 

operations 

possible 

24 

Field 

operations 

not  possible 

0.5 

in. 

rain 

35 

Field 

operations 

not  possible 

26 

Field 

operations 

possible 

27 

Field 

opera tions 

not  possible 

0.3 

in . 

rain 

28 

Field 

opera tions 

not  possible 

0.5 

in. 

ra  in 

29 

Field 

operations 

not  possible 

30 

Field 

opera  tions 

possible 

90- 


October ,  19t6 

• 

Day 

Opera 

.  io  a 

Comments 

1 

Field 

opera tica 

not  possible 

0.2  in.  rain 

2 

Field 

opera l ion 

not  possible 

0.5  in.  rain 

3 

Field 

opera  tica 

not  possible 

0.2  in.  rain 

4 

Field 

operation 

not  possible 

5 

Field 

opera  cion 

net  possible 

6 

Field 

opera  t ion 

po  s  s ib 1 e 

7 

Field 

operation 

not  possible 

0.2  in.  rain 

U 

Field 

opera t ion 

not  possible 

1  in.  snov 

Q 

y 

Field 

opera tioa 

not  possible 

10 

Field 

opera  1 1 0  i  i 

po s 3 lb 1 e 

11 

Field 

opera tion 

possible 

3  no 

12 

Field 

operation 

possible 

13 

Field 

operation 

po s s ib 1 e 

14 

Field 

operation 

possible 

15 

Field 

opera  cion 

possible 

16 

Field 

opera  t ion 

possible 

17 

Field 

opera  tion 

possible 

18 

Field 

operation 

poss ib ie 

19 

Field 

operation 

possible 

20 

Field 

opera tion 

possible 

21 

Field 

operation 

possible 

22 

Field 

operation 

possible 

23 

Field 

operation 

possible 

24 

Field 

operation 

possible 

25 

Field 

opera  tion 

possible 

26 

Field 

operation 

possible 

27 

Field 

operation 

possible 

28 

Field 

operation 

possible 

29 

Field 

opera  tion 

possible 

30 

Field 

operation 

possible 

31 

Field 

operation 

possible 

■ 
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8.3  PROGRAM  USED  TO  ESTIMATE  WORKING  DAYS 
8.3.1  MAIN  PROGRAM 

C  THIS  PROGRAM  IS  USED  IN  CONJUNCTION  WITH  FIVE  SUBROUTINES. 

C  THE  INPUT  DATA  IS  READ  IN  THE  FOLLOWING  FORM 
C  TWO  TABLES  OF  ADJUSTMENT  FACTORS  (Z>  FOR  SOIL  DRYNESS  AS 
C  DESCRIBED  BY  BAIER  AND  ROBERTSON  ARE  READ  ON  CARDS*  20  CARDS 
C  FORMAT  <6X*E5. 2 ) . 

C  MEAN  MONTHLY  EVAPOTRANSP IRATION *  ONE  CARD  PER  STATION  RUT 
C  IF  THE  NUMBER  OF  STATIONS  IS  NOT  SEVEN  THEN  THE  SIZF  OF 
C  PEXBAR  MUST  RE  CHANGED  ACCORDINGLY*  SEVEN  CARDS*  FORMAT (12X* 

C  12F4.2). 

C  READ  COEFFICIENTS  DESCRIBING  UPTAKE  OF  WATER*  8  CARDS  FORMAT 
C  (6F5.2). 

C  READ  CAPACITY  OF  EACH  ZONE*  1  CARD  FORMAT ( 6F5 • 2 ) . 

C  READ  MOISTURE  AT  THE  BEGINNING  OF  EACH  SFASON  AND  CROP 
C  STAGE  DATES*  1  CARD  FORM A T ( 6F 5 . 2  * 9  1 4 ) . 

C  READ  CUMMULATIVE  TOTAL  DAYS  TO  THE  BEGINNING  OF  EACH  MONTH* 

C  1  CARD  FORMAT ( 12 14 ) . 

C  READ  VAPOR  PRESSURES  FOR  FACH  HALF  DEGREE  0  -  110*  111  CARDS 
C  FORMAT (9X*F7.3*3X,F7.5 ) .  THIS  DATA  READ  FROM  TAPE. 

C  READ  SOLAR  RADIATION  AT  THF  TOP  OF  THE  ATMOSPHERE*  365  C  \RDS 
C  PER  DEGREE  OF  LATITUDE.  AGAIN  IN  THIS  IS  NOT  5  THFN  THF  SIZE 
C  OF  OS  MUST  BE  CHANGED..  THIS  DATA  ALSO  RF AD  FROM  TAPE* 

C  FORMAT (72X*F8.1 ) . 

C  READ  X,  DOT*  BLANK*  1  CARD  FORMAT ( 6X * 3A1 ) . 

C  RFAD  MONTHLY  AVERAGE  DAILY  WIND  MILEAGE*  7  CARDS  FORMAT 

C  ( 12F6. 1  ) . 

C  READ  MONTHLY  AVERAGE  HOURS  OF  SUNSHINE*  7  CARDS  FORMAT 
C  (12F6.1). 

C  READ  MONTHLY  AVERAGE  DEWPOINT*  7  CARDS  FORMAT ( 1 2F6 . 1 > . 

C  READ  MONTHLY  AVFRAGE  DA YLENGTH ♦  7  CARDS  FORMAT ( 1 2F6 . 1 > . 

C  READ  TYPE  4  CL IMATALOG I  CAL  CARDS. 

C  OUTPUT  IS  THE  STATION  NO.*  MONTH,  YEAR,  X»S  AND  DOTS  INDICAT- 
C  ING  GOOD  AND  BAD  DAYS. 

DIMENSION  OVW IND(7*1?) *  OVSUN (7*1?) ,  0VDEW(7*1Z>» 

1  A VW I  ND ( 1 2 )  *  AVSUNI12)*  AVDEW(12),  aVDAY(1?)*  QS(5  *365), 
ISM ITTY( 222 ) *  COEF ( 6 )  *MONCON (  12  )  *QCOFF (12) *NUM (11)* 

1 Y ( 6 ) *  V ( 3 1 ) *  KDA  T  F ( 7  )  *  I NFORM ( 80 )  *0 VDA Y ( 7  *  1 2 >  *  Q(365) 
DIMFNSION  PEXBAR ( 7*12) *COEFS( 8*6)  *  CAP AC (6)  *C0NTNT(6  )  * 

1  I  DATES ( 9 ) , DELTA (6 )  *D( 3 ) *  PEME AN ( 12  >  *COF ( 6 ) *  ATARLE ( 100 )  * 

1  BT ABLE (100),  NUMB( 1 0 ) *Z ( 12 ) *CNTENT (6  ) 

INTEGER  P*  DATE*  BLANK*  SNOW 

DATA  LT ,LC  *LL/ ' T  »  * ♦ C  »  *  »L  * / *NUMB/ • 0  *  *  » 1  *  * ♦ 2  *  *  *  3  «  *  »  4 •  * 

1  *  S  •  *  •  6  *  *  •  7  *  *’8*  *  '  Q  *  / 

DATA  KDATE/'APR»*«MAY»*»  JUNE • *  •  JULY •  *  »  AUG*  *  »  SEP T *  *  » OCT  » / 

8  FORMAT ( 6X,5F10.5  ) 

104  FORMAT (I3*I4*3I2*2E3.0*22X»F2.1*A1*1X*F2.1*A1) 

320  FORMAT ( 1X*I4,1X*A4*2H19*I2*2X,31A1) 
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CALL  BEGIN  (  ATABLF,BTABLE,PFXBAR,COEFS*CAPAC,CNTENT  * 

1  I  DATES  *  MON CON  *SMITTY*QS*X*DOT  » BLANK *OVW I ND  *OVSUN *OVDEW * 
10VDAY) 

1010  KNTDAT  =  1 

PRFCIP  =  0.0 
SNOFAL  =  0.0 
DO  97  I  =  1*6 
CONTNT ( I )  =  CNTENT ( I > 

97  CONTINUE 

9  8  READ<4* 104* END  =  501  )  I  ST  A  *  I  STB  *  I YR  ,  I MO  *  I DAY  * TMAX * TM I N > SNO ♦ 

1  SNOW  *  PRF  *  P 
KSTB  =  I  STB 
GO  TO  93 

99  READ(4*l04*FND  =  50l ) ISTA* I STB* IYR  *  I MO  *  I  DAY ♦ TMAX ♦ TMI N  * SNO ♦ 
1  SNOW*  PRE*  P 

93  I  DATE  =  I  MO* 1 00  +  I  DA Y 
ITEST  =  I DATES ( 1 ) 

I E (  I  TEST  .GT.  IDATE)  GO  TO  99 
I E ( I STB  .NE.  KSTR)  GO  TO  1010 
KNTDAT  =  1 
DO  50  I  =1*31 
V (  I  )  =  BLANK 
50  CONTINUE 
MN  =  1 

CALL  START (  I  STB  * AVW I ND  *OVW I ND  *  AVSUN  *OVSUN  *  AVDEW  *OVDEW  * 
1AVDAY*0VDAY*PEMEAN  *PEXBAR *Q *QS ) 

M  =  0 
KMO  =  I  MO 
GO  TO  130 

100  READ(4*104*FND  =  5  01 ) ISTA*ISTB* I YR *  I  MO  *  I  DA Y * TM AX  * TM I N ♦ SNO * 
1 SNOW  *  PRE,  P 

130  IDATE  =  I  MO* 100  +  IDAY 
M  =  M  +  1 

I F ( I MO  .NE.  KMO)  GO  TO  302 

131  I F ( I  MO  .GE.  11)  GO  TO  600 

101  IF (  ITEST  .GT.  IDATE)  GO  TO  81 
DO  80  I  =  1*6 

COEE(I)  =  COEFS( KNTDAT, I ) 

80  CONTINUE 

KNTDAT  =  KNTDAT  +  1 
ITEST  =  I  DATES ( KNTDAT ) 

81  SUMDEL  =  0.0 
KMO  =  I  MO 

AMDEW  =  AVDEW (  I  MO ) 

IF (P.EQ.LT.OR.P.EQ.LC.OR.P.EQ.LL )G0  TO  1110 
DO  17  I  =  1*10 
I E  (  P  .EQ.  NUMB ( I ) )  GO  TO  18 

17  CONTINUE 

I F  ( P  .EQ.  BLANK)  GO  TO  1110 

18  PRECIP  =  PRE  +  (I  -  1)*  0.01 
1110  CONTINUE 
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CALL  E VAPOR ( TMAX*TMIN*SMITTY*  AMDEW  *  MON CON ♦ I  MO  *  I  DAY  >Q* 

iavsun*avday,avwind*el  ) 

CALL  MOIST ( PEMFAN  *CONTNT  *CAPAC*Cl *C2  *COEF  *KNTDAT *NEW* 
1ATABLE*BTABLE *PRFCIP *D*EL* IMO  ) 

PRECIP  =  0.0 
SNOFAL  =  0*0 

IF ( SNOW.EQ.LT. OR. SNOW. EQ.LC. OR. SNOW. EQ.LL)  GO  TO  20 
I F ( SNOW  .FQ.  BLANK)  GO  TO  20 
DO  21  I  =  1*10 

IF (SNOW  .EQ.  NUMB ( I ) )  GO  TO  22 

21  CONTINUE 

22  SNOFAL  =  SNO  +  (I  -  1)*0.01 
20  CONTINUE 

IF ( I MO  .FQ.  4  .AND.  TMAX  .LT.  45.0)  GO  TO  301 
I F ( SNOFAL  .GT.  0.0  )  GO  TO  304 

I F ( I  MO  .LE.  5  .AND.  TMAX  .GT.  45.0)  SNOWFL  =  0.0 
I F ( I MO  .GE.  6  .AMD.  TMAX  .GT.  40.0)  SNOWFL  =  0.0 
I F ( SNOWFL  .GT.  0.0)  GO  TO  301 

IF (CONTNT ( 1 ) .GE. 0.190. OR. CONTNT ( 2 ) .GF. 0.285 )  GO  TO  301 
V  (  M  )  =  X 
GO  TO  100 

301  V ( M )  =  DOT 
GO  TO  100 

304  SNOWFL  =  SNOFAL 
GO  TO  301 

302  DATE  =  KDATE(MN) 

WR I TF ( 8  *  320 )  ISTR,  DATE  ♦  I YR  *  V 
MN  =  MN  +  1 
DO  303  I  =  1*31 
V (  I  )  =  BLANK 

303  CONTINUE 
M  =  1 

KMO  =  I  MO 
KSTB  =  I  STB 
GO  TO  131 
600  CONTINUE 

CALL  WINTER (CONTNT  *CAPAC ♦ I YR , I  STB ) 

GO  TO  99 
501  CONTINUF 
REWIND  4 
END  FILE  8 
REWIND  8 
STOP 
END 
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8.3.2  SUBROUTINE  TO  READ  INPUT  DATA 


2 

3 

4 

9 

12 

14 

16 

91 

107 


45 


13 


30 

19 


20 


21 


22 


23 


SUBROUTINE  BEGI N ( AT ABLE  * ST ABLE  * PEXBAR * COEFS * CAP AC  * CNTEN T * 
1  I  DATES  *  MON CON  * SMI  TTY *QS  *X  *DOT  * BLANK  *OVW  I  ND  *OVSUN *OVDEW  ♦ 
10VDAY ) 

DIMENSION  AT  ABLE ( 100)  *  BT  ABLE ( 1 00 )  * P FXRAR ( 7 * 1 2 > *C0EFS(8*6> 
1  *  CAP AC ( 6  )  *  CNTENT ( 6 )  *  I  DATES (9)  *MONCON( 12 )  * SM I TT Y ( 2 2 2  >  * 

IQS (5  *365) *OVW I ND ( 7  * 1 2 )  *OVSUN(7*12)  *OVDFW(7*12) *OVD AY 
1(7*12) 

FORMAT ( 6F5.2  *9 14  ) 

FORMAT ( 12X,12F4.2> 

FORMAT ( 1X*6F10.2 ) 


FORMAT (6X*10F5.2) 

FORMAT ( 9X*F7.5*3X*F7.5) 

FORMAT (72X*F8.1  ) 

FORMAT ( 1 2F6 • 1 ) 

FORMAT (6X*1A1*A1 *A1 *  ) 

FORMAT (1214) 

RE AD ( 5  *  9 )  ATABLE,  BT  ABLE 
1*7 

( PEXBAR ( I  * J ) *  J=l*12) 


DO  45  I  = 
READ ( 5  *  3 ) 
CONTINUE 
DO  13  I  = 
READ ( 5  *  2 ) 
CONTINUE 
READ ( 5*2) 
RF  AD ( 5  *  2 ) 


1*8 

( COFFS (  I  * J) 


J  =  1*6  ) 


CAPAC 

CNTENT*  IDATES 
READ (5*107)  (MONCON(I)  *  1=  1*12) 
RF  AD ( 3  *  1 2 )  ( SMI  TTY ( I )  *1=1  *222  ) 

DO  19  I  =  1*  5 
DO  30  J  =  1*365 
READ(3*14)  QS ( I  *  J ) 

CONTINUF 
CONTINUE 
REWIND  3 
READ ( 5  *  9 1  ) 

DO  20  I  = 

RE AD ( 5  *  1 6 ) 

CONTINUE 


X*  DOT,  BLANK 
1*7 

( OVW I ND ( I  * J) * 


J=1 *12} 


WR I T  F ( 6  *  4 )  OVW I ND (7*12) 

DO  21  I  =  1*7 

RE  AD ( 5  *  1 6 )  ( OVSUN ( I  * J  )  *  J=l*12) 

CONTINUE 

WRITE (6 *4)  OVSUN (7* 12) 

DO  22  I  =  1*7 

RE AD ( 5  *  1 6 )  ( OVDEW ( I  ,  J )  ,  J= 1*12) 

CONTINUE 

WR  I  TE ( 6  *  4 )  OVDEW ( 7  *  1 2 ) 

DO  23  I  =  1*7 

RE AD ( 5  *  1 6 )  ( 0 VDAY ( I  *  J ) *  J=l*12  ) 
CONTINUE 

WR  I  T E ( 6  *  4 )  OVDAY (7*12) 

RETURN 

END 
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8*3.5  SUBROUTINE  TO  SET  INITIAL  CONDITIONS 

SUBROUTINE  START ( I  STB  *  AVWI ND  >OVW IND  * AVSUN  *OVSUN  *  AVDEW  * 
10VDEW*AVDAY,0VDAY*PEMEAN ,PEXBAR*Q*QS) 

DIMENSION  AVWIND(12)  *  OVW I ND ( 7  *  1 2 )  *AVSUN(12) *OVSUN(7*12)  * 
1 A VDEW (12) ♦OVDEW(7*12) *AVDAY( 12 ) *OVDAY(7*12> *PEME AN ( 12 > * 
1PEXBAR( 7*12 ) *Q(365  )  *  QS(5*365) 

I E ( I  STB  *EQ.  4480  )  1  =  1 

I E (  I  STB. EQ. 3 88  0. OR. I STB.EQ. 3890 >  1=2 

IE ( I STB.EQ .1080. OR .ISTB.EQ.1085.OR.IST8.EQ.il 00 )  I  =3 


IF ( ISTB 

.  EQ  . 

2195 

•OR.  ISTB. 

EQ.  2200)  I  =  4 

IE (  ISTB 

.EQ. 

2205 

.OR.  ISTB 

•EQ.  2208)  I  =  4 

I E ( ISTB 

•  EQ. 

5168 

.OR.  ISTB 

•EQ.  1820)  I  =  4 

IF ( ISTB 

•  EQ. 

6800 

)  I  =  5 

IE ( ISTB 

.EQ. 

2520 

)  I  =  6 

IE ( ISTB 

•  EQ. 

0560 

)  I  =  7 

DO  41  J 

=  1* 

12 

AVWIND(J)  = 

OVW I ND ( I  * J) 

CONT INUE 
DO  42  J 

=  1* 

12 

AVSUN ( J )  =  OVSUN ( I  * J ) 

42  CONTINUE 

DO  43  J  =  1*  12 
AVDEW ( J )  =  OVDEW ( I , j) 

43  CONTINUE 

DO  44  J  =  1*  12 
AVDAY ( J )  =  OVDAY (  I  * j) 

44  CONTINUE 

DO  46  J  =  1*12 
PEMEAN(J)  =  PEXBAR (  I  *  J  ) 

46  CONTINUE 

IE ( I  STB. EQ. 44  80. OR. I STB.EQ. 3880 .OR . I  STB. EQ • 3890 >  1  =  1 

I E ( ISTe.EQ.1080.0R. I  STB . EQ . 1 085 .OR .  I  STB . EQ . 1 1 00 )  I  =  2 
IF( I STB.EQ. 2 200. OR. I  STB . EQ . 2 195 .OR . I  STB • EQ. 6800 )  1=3 

IE ( I  STB • EQ. 2205  .OR.  ISTB  .EQ.  2208>  I  =  3 
I F ( I STB  .EQ.  5168  .OR.  ISTB  .EQ.  1820)  I  =  3 
IF ( I STB  .EQ.  0560  )  I  =  4 
I E (  ISTB  .EQ.  2520  )  I  =  5 
DO  40  J  =  1*365 
Q( J)  =  QS ( I  * J ) 

40  CONTINUE 
RETURN 
END 


-  p  p  ~ 


.  T  T  r  n-im  '  •  T  T  T  VO  T  '  ~  T  •  *  T  t  q  :  -■  c*  .  .  ' 

•  Pl*IWVA*«TP 

'  q  .  0 . 5  '  >  Qe  '  '■  •  v  -,  j  ,  1 A  .  A  V  A  ,  •  ->r»VC  I 
•  •  .  .  -  s .  ;  f  •  ■  )•  ( Sf  JOHI  WVA  70 !  ' 

(  '  p  •  j )  p  o  .  (  <=»  <>  p )  .  ( 5  r  •  •  )  ?  \  p  *  -i  q  r 


r  = 

1  ( 

?A  AA 

.  . 

3T?T ) 3: 

5=1 

(A'Apf.  3. 

q  T  r  T  . 

.PQ  .  " 

.pTr n q ’ 

.  1  P  f  . 

•  . 

-  D  r  . 

. 

A  =  !  (^c<; 

,  V’.  Tq  I 

•  AO. 

ap  f  <; 

• .  • 

^ TP  I  )  H 

A  =  I  lan^ 

.O'3.  r>T?T 

.*  . 

*P55 

•  0  3  , 

,  T  T  )  i  T 

A  =  T  ( "5*  r 

•O'.  QT'I 

. 

M  r  * 

. 

ATP  T ) 3! 

? 

*  T  ( 

"PA  3 

.  >  . 

PT  -»T  )  ;r 

P 

*  I  ( 

C5£5 

.03. 

8  T  P  T  )  *1 

X 

«  I  ( 

oppp 

.  3. 

•'TP  T  )  ■  T 

5  r 

* 

l  r  a  io 

<  L. T ) 0^  r ' VO  =  < LIT  M  VA 

‘-UWITtfOO  r  A 

51*  f  •  L  5A  00 

(L.pnipvo  =  (ur-pv". 

3UKITK0*)  sa 
Sf  *r  *  L  PA  00 
<L-  T  )  ’«n  '"'■  ■')=  (  L  )  \H  p*  V A 

^  3lWlTH0:>  PA 

51  *  I  ~  L  A  A  0( 

(  L  •  I  )  Y  AG  VO  *  (  L  )  VA  '’VA 

A  A 

5  r  •  r  =■  L  PA  or 
<L*  I  )  "I A  3*  *  (L)HA  ^'/^q 

lU^ITKIOD  PA 

•TC.  T^.  'A/Oflf,  q  •  q  T  ■*  I  •  .O-  T  )  ’  I 

.  f  .  .  ’  .  .  .  '  HI 

.  -  '  ■ .  .  r 5#  ' .  *  ■  ■ .  ■  »oo5  .  '  . c: ^ n q t 

r  =  I  V'^CS  .  .  3  TP  I  .qr. 
p  =  i  ('>5*r  .03*  ct q t  .or.  «pt*  *03#  flTp.nii 

A  «  T  (  Od?n  ,0'U  «T’T)3T 
*  *  T  <  P535  #0  .  pTP T ) 31 
PPP.f  =  L  0 A  OG 

<l*  neo  =  (do 

^Ui*IT*Q“>  r  A 

«P!jT  3* 


-96- 


8.3.4  SUBROUTINE  TO  FIND  SOIL  MOISTURE 

SUBROUTINE  MO  I ST ( PEME AN * CONTNT * CAP AC  * C 1  *  C2  *COEF  *  KNTDAT * 
1NEW*ATABLE*RTABLE *PREC I P *D *EL *  I  MO  ) 

DIMENSION  Y ( 6 )  *COF(6)  *  DELTA (6) * PFMF AN ( 1 2 > * CON TN T ( 6 )  * 

1  CAP AC (6 )  *COEF(6  > * ATARLF < 100  >  *BT ABLE ( 100 ) *D( 3 } 

PE  =  EL  *  .0034 
RUNOFF  =  0.0 
DO  15  I  =  1*6 
Y ( I )  =  0.0 
15  CONTINUE 

PEDIF  =  PE  -  PEMEAN ( I  MO ) 

AYE  =  0.0 

DO  110  I  =  1*6 

DEL  =  CONTNT (  I ) /CAPAC (  I  ) 

W  =  7.81  -  0. 11*DEL*100. 

I F ( W  .GT.  0.)  GO  TO  162 
W  =  0.0 

162  COF (  I  )  =  COEF (  I  ) 

I F ( I  .LE.  3)  GO  TO  98 
DO  96  J  =  2*1 
K  =  J  -  1 

COF ( I )  a  COF ( I )  +  COF ( I )*COF(K)*( 1.0  -  CONTNT ( K) /CAP AC 
1  (K)  ) 

96  CONTINUE 
98  IT  =  DEL*100. 

IF (IT  .GT.  0)  GO  TO  102 
WORK  =  0. 

W  =  0  • 

GO  TO  108 

102  I F ( KNTDAT  .GT.  D  GO  TO  106 
WORK  =  ATABLE ( IT ) 

GO  TO  108 

106  WORK  =  BTABLE ( IT  ) 

C  FIND  ACTUAL  FVAPOTRANSP I R AT  I ON 

108  IT  =  ( COF ( I ) *DEL*WORK*PE*EXP (-W*PED I F  )  +  .005)*100. 

DEL  =  IT 

DEL  =  DEL*. 01 

I F ( DEL  .LE.  CONTNT ( I ) ) GO  TO  109 
DEL  =  CONTNT ( I ) 

109  DELT A ( I  )  =  DEL 

110  CONTINUE 

IFCPRECIP  .EQ.  0.)  GO  TO  220 
IF ( PREC I P  .GT.  1.0  )  GO  TO  128 
AYE  =  PREC  I P 
GO  TO  140 

C  DETERMINE  INFILTRATION 
128  OGER  =  ALOG ( PREC I P ) 

AYE  =  0  #9177+  1.811 *OGER  -  0 . 0097*OGER*CONTNT ( 1 ) /CAPAC ( 1 ) 
IT  =  (AYE  +  . 005 ) *  1 00 • 

AYE  =  IT 
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AYE  =  AYE*. 01 
RUNOFF  =  PRECIP  -  AYE 
I F ( RUNOFF  .GE.  0.)  GO  TO  150 
140  RUNOFF  =  0. 

C  DETERMINE  MOISTURE  DEFICIENCY 
150  DO  180  I  =  1*6 

DEFCIT  =  CAP  AC ( I )  -  CONTNT(I)  +  DELTA (  I ) 
IF (AYE  .GT.  DEFCIT)  GO  TO  170 
Y  (  I  )  =  AYE 
AYE  =  0. 

GO  TO  220 
170  Y ( I )  =  DEFCIT 

AYE  =  AYE  -  DEFCIT 
180  CONTINUE 

C  DETERMINE  CHANGE  IN  MOISTURE  CONTENT 
220  DO  230  I  =  1*6 

CONTNT ( I )  =  CONTNT ( I )  -  DELTA(I)  +  Y ( I ) 
IF ( CONTNT ( I )  .GE.  0.0)  GO  TO  230 
CONTNT ( I )  =  0. 

230  CONTINUE 
RETURN 
END 


8.3.5  SUBROUTINE  TO  FIND  EVAPORATION 

SUBROUTINE  EVAPOR ( TM AX , T M I N  * SM I  TTY  * AMDEW *MONCON *  I  MO* 

1 IDAY*Q*AVSUN*AVDAY*AVWIND*EL) 

DIMFNSION  SM I  TTY ( 222  *  *MONCON( 12 >  *0(365 )  *AVSUN(12) * 
1AVDAY(12)>  AVWIND(12) 

1110  I F ( TMAX  .LT.  32.0  )  GO  TO  998 

IF (TMAX*TMIN  .LE.  0.)  GO  TO  998 
INDEX  =  TMAX  +  TM I N  +  1.5 
AIR  =  SM I  TTY ( INDEX ) 

INDEX  =  2.0  *AMDEW  +  1.5 

AIR  =  (AIR  -  SMITTY( INDEX) )*33. 86389 

INDEX  =  MONCON ( I  MO )  +  IDAY 

QO  =  Q ( INDEX) 

SUN  =  AVSUN(IMO) 

DAY  =  AVDAY ( IMO) 

WIND  =  A VW I ND ( IMO) 

QSS=  QO* (0.251  +  0 • 6 1 6* ( SUN /DAY ) ) 

EL  =  -53.39  +0  •  337*TMAX+  0 . 5 3  1  * ( TMAX-TM  I N )  +  0.010  *QO  + 
10.051 2*QSS  +  0.09 77* WIND  +  1.77*AIR 
I F ( EL  .GT.  0.0)  GO  TO  37 
998  EL  =  0.0 
37  CONTINUE 
RETURN 
END 
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8.3.6  SUBROUTINE  TO  DETERMINF  AMOUNT  OF  WINTER  PRECIPITATION 
CONSERVED 


104 


601 


602 

603 

604 


605 


612 


613 

615 

620 


608 

607 

609 


610 

501 


SUBROUTINE  W I NTER ( CON  TNT  * CAP AC  * IYR.ISTB) 

INTEGER  P 

DIMENSION  CON  TNT (6 )  *NUMB (10)  *  CAP AC ( 6 )  *Y (6 ) 

FORMAT ( 13*14*312*2 F3.0*26X*F2.1*A1) 

DATA  LT*LC*LL/'T»*»C«*'L'/ *NUMB/»  0 **,1»*‘2'*‘3‘**4»* 
l,5,*,6,*,7,*,8,*,9,/ 

PRECIP  =  0.0 

READ (4*1 04*END=501  )  I  ST A  *  I  STB  *  I YR *  I  MO  *  I  DAY  * TM AX  * TM I N * 
1 PRF  *  P 

IDATE  =  I MO*l 00  +  IDAY 

IF (P.EQ.LT •OR.P.EQ.LC.OR.P.EQ.LL)  GO  TO  604 

DO  602  I  =  1*10 

IF ( P • EQ •  NUMB ( I ) )  GO  TO  603 

CONT INUE 

PRE  =  PRE  +  (  I  -  1 >*0.01 
PRECIP  =  PRECIP  +  PRE 
I E ( I  DATE  .EQ.  331)  GO  TO  605 
GO  TO  601 

I F ( I  ST  A  .GE.  303  )  GO  TO  612 
AYE  =  PREC I P*0 .25 
GO  TO  615 

I E ( I STA  .EQ.  307)  GO  TO  613 

AYE  =  PREC IP*0.35 

GO  TO  615 

AYE  =  PREC I P*0  • 1 0 

DO  620  I  =  1*6 

Y( I)  =  0.0 

CONTINUE 

DO  607  I  =  1*6 

DEECIT  =  CAP AC ( I )  -  CONTNT(I) 

I E ( AYE  .GT.  DEFCIT)  GO  TO  608 

Y ( I )  =  AYE 

AYE  =  0.0 

GO  TO  609 

Y ( I )  =  DEFCIT 

AYE  =  AYE  -  DEFCIT 

CONTINUE 

DO  610  I  =  1*6 

CONTNT (  I  )  =  CONTNT (  I  )  +  Y(I> 

I F ( CONTNT ( I )  .GE.  0.0)  GO  TO  610 

CONTNT (  I  )  =  0.0 

CONTINUE 

CONTINUE 

RETURN 

END 
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